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Abstract 

The properties of the quark-gluon medium observed in high energy 
nucleus-nucleus collisions are discussed. The main experimental facts 
about these collisions are briefly described and compared with data 
about proton-proton collisions. Both microscopic and macroscopic 
approaches to their description are reviewed. The chromo dynamics of 
the quark-gluon medium at high energies is mainly considered. The 
energy loss of partons moving in this medium is treated. The principal 
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conclusion is that the medium possesses some collective properties 
which are crucial for understanding the experimental observations. 

A foreword 

It is a great honor to us to publish this paper in the volume dedicated 
to the memory of V.L. Ginzburg. His contribution to the condensed matter 
physics crowned by the Nobel Prize is extremely important. Nowadays these 
ideas become actual for high energy physics. The properties of the quark-gluon 
medium created in high energy nucleus collisions ( often called the quark-gluon 
plasma) are intensively studied. The first steps of chromodynamics of quark- 
gluon matter strongly remind those of electrodynamics of continuous media 
with some visible differences, however. 



1 Introduction 

The topic of the quark-gluon medium is so widely spread nowadays that it 
deserves a series of monographs, even in spite of still being in the actively 
developing stage, and can not be fully covered within a single review paper. 
That is why we mostly concentrate here on general ideas about the evolution 
of the quark-gluon medium in high energy heavy ion collisions and on those 
properties of the medium which are revealed by the energy losses of partons 
moving in it. They are described by the chromodynamics of the quark-gluon 
medium which is of the main topic of this review. 

Nevertheless, it would be too "ostrich-like" behavior not to mention the 
nearby problems and approaches. Therefore we include in the present re- 
view brief discussions of some of them. In particular, the analogies and 
differences between electrodynamical and chromodynamical processes are de- 
scribed. The extremely popular hydrodynamical approach has been widely 
discussed in many review papers. That is why we decided to give just a short 
description of it here referring to the corresponding literature. The Anti de- 
Sitter/Quantum Chromodynamics (AdS/QCD) correspondence is touched 
upon only slightly. Lattice results are mentioned very briefly too. They are 
mostly applicable to static properties of hadrons and are of interest to us 
here in what concerns phase transitions. The very interesting ideas about 
possible local CP violation in heavy ion collisions are briefly mentioned but 
not discussed in detail. The relation of the little bang in heavy ion collisions 
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to the Big Bang in cosmology is another topic of interest. However, it is out 
of the scope of our survey. A short overview of experimental results is pre- 
sented in section 2. However, in subsequent sections of this review only those 
experimental papers are cited which are discussed in more details in connec- 
tion with theoretical considerations. We apologize to those whose work could 
not be completely reviewed due to the limited size of this paper. 

The internal structure and interactions of particles and nuclei are the 
main goal of studies at accelerators. According to present knowledge strong 
interactions are described by quantum chromodynamics (QCD). It describes 
quarks and gluons as elementary objects (partons) responsible for the inter- 
action. 

At low energies the partons are confined inside strongly interacting par- 
ticles (hadrons) and define their static properties. At large densities we have 
the non-perturbative region of QCD describing the strongly interacting mat- 
ter in thermal equilibrium at finite temperature. 

The matter produced in collisions is surely different from the above- 
described one. Highly coherent parton configurations and strong internal 
fields become especially important. By colliding two heavy nuclei at ultra- 
relativistic energies one expects to get a hot and dense internally colored 
medium. It should exhibit some collective properties different from those 
seen at static conditions. Low-momenta modes in nuclear wave functions 
can be described in terms of classical fields coupled to some high-momenta 
static color sources. 

According to present knowledge, the collision process (energy evolution 
of all modes) proceeds through several stages. We discuss some of them at 
the beginning of the review. The corresponding equations describing these 
processes are also discussed. At one of the stages quarks and gluons may 
become deconfined inside some finite volume during a short time interval. 
According to present knowledge they form some ideal liquid. Its collective 
properties are revealed by its mechanical motion as a whole described by 
hydrodynamics and by chromodynamical response to partons penetrating it 
described by in-medium QCD. 

At large transferred momenta the QCD coupling becomes smaller (the 
asymptotical freedom) and the perturbative approach starts being applicable. 
Experiments on high energy particle collisions are needed precisely to make it 
possible to study the processes with high transferred momenta. In accordance 
with the uncertainty principle it implies a possibility to learn the particle 
structure at ever smaller distances. In this way the quark-gluon content of 
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a particle, its scale dependence and the properties of the interaction region 
are studied. 

Moreover, the hard process inside medium may induce its collective co- 
herent response. That would allow to study medium properties like explo- 
sions in earth depths serve to learn its structure. This is especially true at 
highest energies because the parton densities inside colliding strongly inter- 
acting nuclei increase with energy. Therefore, the collective effects become 
pronounced and, in central heavy ion collisions, the long-range correlations 
appear alongside with the short-range correlations . 

To learn experimentally about the properties of the matter inside the in- 
teraction region one should study the energy losses of various trial partons 
moving in it. This is analogous to studies of energy losses of electrons pass- 
ing through the ordinary amorphous matter. Recent experimental results of 
SPS (Super Proton Synchrotron) and RHIC (Relativistic Heavy Ion Collider) 
clearly show that the collective medium behavior becomes crucial in heavy 
ion collisions. The radiation properties of partons are modified and energy 
spectra of final hadrons as well as their correlations are changed relative to 
those in pp collisions. The induced coherent emission reveals the collective 
response of the medium to partons penetrating it. The medium itself dis- 
plays the collective motion. The hadrochemical content of produced particles 
changes. Both microscopic and macroscopic QCD approaches to theoretical 
description of these findings are effective. Mechanical and thermodynamical 
properties of the medium are studied in the hydrodynamical approach. We 
describe all of them in what follows together with discussion of corresponding 
experimental data. Everywhere in this review we use the system of units in 
which h — c—1. 

2 The main experimental findings 

The nucleus- nucleus (AA) collisions are studied experimentally at SPS and 
RHIC. The very first question addressed is the difference in AA and pp 
processes at high energies. Various characteristics were measured and com- 
pared. They clearly show that nuclear collisions can not be considered as 
an incoherent superposition of the nucleon-nucleon collisions. Some collec- 
tive properties of the medium must be taken into account to explain the 
difference. 
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The yield of particles produced at rapidity y and large transverse mo- 
menta pt is suppressed in AA relative to pp. This is demonstrated by plotting 
the nuclear modification factor Raa for single-inclusive hadron cross section: 

c / 7 \ dN AA /dydp T 

RAA{y,pr,b) = — — j- ■ 1 
NpartdNpp/dydpT 

The factor Raa measures the deviation of AA reaction with the number 
of participant nucleons N part (determined by the nuclear overlap function 
at impact parameter b) from an incoherent superposition of pp reactions 
for which Raa = 1- A similar expression can be written for the dihadron 
correlations (usually called Iaa)- 

Both in particle and heavy ion collisions the production of the soft mo- 
mentum particles with the transverse momenta pt < 2 GeV strongly domi- 
nates. The corresponding spectrum decreases exponentially with transverse 
kinetic energy y / pfT+ mf — rrii. The pion spectra are flatter in the AA colli- 
sions than in the pp ones. This is even more true for kaons and antiprotons 
with the latter having the smallest slope. These spectra are used to deter- 
mine the kinetic freeze-out temperature and the dependence of the transverse 
flow velocity on the energy and centrality (i.e. degree of collectivity of the 
interaction). 

Accordingly, soft pions dominate multiplicity distributions. Their shape 
is approximately described by the negative binomial distributions both in 
particle and heavy ion collisions at high energies albeit some deflections from 
them (predicted by QCD) are revealed by such correlation measures as their 
moments. 

The shape of the charged hadron transverse momentum spectra changes 
at p T ps 3 GeV from the exponential to the power-like as predicted by the per- 
turbative QCD. This power law tail affects very small yield of particles. The 
RHIC results on single particle inclusive distributions in central Au-Au colli- 
sions at the energy of 200 GeV show the strong pr-independent suppression, 
Raa ~ 0.2, of the hadronic yield at large px > 4 GeV. The measurements 
covered the interval in p T up to 20 GeV. This striking deficit of the high 
transverse momentum particles points out to the energy loss of partons in 
the medium. It corresponds to the so-called jet quenching effect with soft- 
ened spectra of hadrons produced by partons in the medium compared to 
the vacuum. Therefore, this suppression factor is a powerful tool to map out 
the density of the medium. 
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The high p T suppression in Au-Au and Cu-Cu collisions is comparable at 
the same number of "participants" N part (nucleons which underwent at least 
one inelastic interaction). The ratio Raa is closer to 1 for more peripheral 
collisions. It depends on the angular orientation of the high-p^ hadron with 
respect to the reaction plane (on the azimuthal angle). This dependence 
may be explained as due to the dependence of the energy loss on the length 
traversed by the particle in the medium. 

The ratio Raa does not depend on hadron identity at high pr- In contrast 
to this, in electrodynamics the radiation of electrons is much stronger than 
that of the muons. Interesting effects take place for suppression of the spectra 
of charmed mesons and charmonium. The former are suppressed at the same 
level as pions while the latter is suppressed at the same level as at the SPS 
energy y/s = 17 GeV. 

Another intriguing feature of the data is the fact that direct photons seem 
to be not suppressed in the region of p T from 4 GeV to 15 GeV while at very 
large ~ 20 GeV preliminary data shows a suppression with Raa ~ 0.6. 

At smaller transverse momenta, 0.25 < px < 4 GeV, the suppression 
effect is weaker, Raa even increases slightly at low px < 3 GeV. This is qual- 
itatively explained by the parton rescattering in the initial state which leads 
to the widening of the transverse momenta spectra (Cronin effect). At the 
SPS energies the results refer only to this domain of transverse momenta with 
larger values of Raa than at RHIC. The SPS value (Pb-Pb at ^=17.3 GeV) 
of Raa at = 4 GeV is much larger than at higher energies of RHIC exceed- 
ing 1 for pt > 2 GeV (strong Cronin effect). That is awaited from the more 
important role of gluons as compared to quarks at RHIC (C A /C V = 9/4). For 
collisions with smaller overlap d-Au at 200 GeV the Cronin effect persists at 
least up to pt=8 GeV. Experimental data show that the suppression of inclu- 
sive spectra sets in somewhere in between ~ 20 GeV and \/s ~ 60 GeV. 
The dominating process of soft particle production is characterized by the 
scaling with the number of participants rather than by that in the number 
of binary collisions typical for hard processes. 

Correlations among the secondary particles are the more sensitive in- 
strument in distinguishing between models predictions than the single par- 
ticle inclusive data. Their studies revealed many exciting features. Most 
spectacular of them refer to jets produced in the processes in which high 
Pt particle is present. If a high-p T particle is chosen as a trigger and 
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other hadrons are measured within a surrounding it circle^] defined by the 
(pseudo)rapidity+azimuthal angles in the plane orthogonal to the momen- 
tum of the initial particle, then the Gaussian-like structure around the trig- 
ger (the near-side jet) is usually seen in pp-reactions. Moreover, both two- 
and three-particle correlations show two back-to-back jet-like peaks (di-jet 
production). Theoretically, jets are considered, at a first approximation, as 
survived remnants of hard-scattered quarks and gluons. 

In AA collisions, jets are expected to widen in angles and momentum 
spectra of their particles are expected to soften. This phenomenon is termed 
jet quenching. Quantitative analysis of jet quenching in heavy ion collisions 
requires model building. So far there are many uncertainties in their pre- 
dictions. Also, the biases are introduced by some experimental jet selection 
criteria. To be less dependent on these effects, different triggers (e.g., 7 and 
7T°) were used and the away-side hadron spectra were measured. 

Correlations of non-photonic electrons from charm and bottom decays 
with hadrons are different for jets with D and B mesons. 

Multiparticle correlations and fully (calorimetrically) reconstructed jets 
are the main goals of recent efforts in the studies of pp and AA collisions. 
First results for fully reconstructed jets in pp, Cu-Cu, Au-Au collisions clearly 
demonstrate that jets are broadened in the quark-gluon medium. In general, 
The color connections of jets to the bulk or, e.g., to beam remnants in W+jet 
events and the event-by-event analysis of color flows are also important. The 
first experimental hints of intra-jet broadening due to such color flows in 
them are obtained. 

The influence of jet quenching and background fluctuations on conclu- 
sions about jet fragmentation functions is actively studied now. Surprisingly 
enough, no modification of them in AA collisions compared to the pp ones 
was observed at high in preliminary studies. Probably, it shows that a 
biased jet population was selected (surface effect?) and the criteria of jet 
selection should be modified. 

Even more surprising were observations of completely new topologies 
called ridge and double-humped events in AA collisions. In central colli- 
sions the near-side jet happens to stand on the pedestal (ridge) with long 
extension in measured pseudorapidities up to \r]\ =4 and steeply falling in 

3 Beside such cone algorithm, the recombination algorithms with iterative pairing of 
nearby particles and simple Gaussian filtering are very widely used. These algorithms are 
helpful in the separation of jets from the background. 
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the azimuthal direction. The presence of the ridge is independent of the jet 
peak. The characteristics of this peak are the same as those of bulk particle 
production. However, the particle spectrum in the ridge is somewhat harder 
than that in the bulk. Moreover, this structure is also seen in two-particle 
correlations without a jet trigger. These events include all charged hadrons 
with p T > 150 MeV. The event-by-event three-particle correlation studies in 
central collisions at higher > 1 GeV showed that the particles from the 
ridge are uncorrelated in pseudorapidity between themselves and with parti- 
cles in the peak. The production of the ridge is uncorrelated with the jet-like 
particle production. The ridge disappears in the low-multiplicity peripheral 
collisions and at high p T of the trigger particle. Both the long pseudora- 
pidity extension of the ridge and large wide clusters found in two-particle 
correlations points out to the important role of the collective effects. 

The away-side peak at A0 = 7r is present in pp-collisions but is replaced 
by the broad away-side structure in Au-Au collisions. The two symmetrical 
maxima (humps) are clearly seen in most central collisions at A0 it ± 1.1. 
They are resolved after the flow effects are subtracted and contain hadrons 
with small energies. The peak positions are approximately independent of 
the transverse momenta of the trigger and associated particles. The peak 
heights are slowly increasing with p 1 ^ 9 at fixed p^ ssoc . The two humps seem 
to merge in a single wide hump at large p T of the trigger (6 < p% %9 < 10 GeV). 
It indicates on the reappearance of the away-side jet just in-between these 
humps as it should be in the finite size medium where the parton with extra 
Pt goes outside it in the form of a jet. The heavy quark jets possess the similar 
qualitative features, albeit measured with the lower statistics. These peculiar 
features are seen both in two- and three-particle correlations. They evidence 
for a conical emission pattern. The existence of these features is surely related 
to the collective properties of the medium. In mid-central collisions the 
maximum position happens to be shifted for some special choices of the 
trigger location relative to in- and out-of- planes directions. This special 
feature can be explained by the effect of the wake (see below). 

The associated yield of away-side hadrons is suppressed at large pt in 
Au-Au collisions relative to the pp ones, albeit less than the single particle 
spectrum (Iaa ~ 0.35 — 0.5 while Raa ~ 0.2). 

Two-particle correlations between charged hadrons generated by final 
state Coulomb interactions provide valuable information. Also Bose-Einstein 
(BE) correlations between pairs of identical bosons caused by quantum sym- 
metrization of wave functions are studied with the help of a'la Hanbury 
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Brown-Twiss (HBT) interferometry. They are observed as an enhancement of 
pairs of same-sign charged pions at small relative momenta and tell us about 
the space-time structure and evolution of the source. In central collisions at 
midrapidity the three radii Ri ong , Rout, Rside are usually reconstructed from 
measurements. They are directed, correspondingly, along the beam direc- 
tion, that of the emission of the pair p T and the one perpendicular to these 
two. Within experimental errors all the three radii turn out to be equal 
and rather steeply decreasing with transverse momenta. This contradicts 
the hydrodynamical predictions and is, therefore, called the "RHIC HBT 
puzzle" . Moreover, recent data show that the HBT characteristics of pp col- 
lisions behave in the similar way to those of AA which adds new questions 
to this puzzle. The size of the correlated particle emission region increases 
with the particle multiplicity in the event. The number of the radii is larger 
(6) for non-central collisions. They grow with the number of participants 
and exhibit oscillations as functions of the azimuthal angle. The growth rate 
decreases with transverse momenta (as well as the radii themselves). The 
amplitude of oscillations is larger for more peripheral collisions and vanishes 
in central collisions. 

The modified shape of bosonic resonances created in nuclear collisions 
has been observed. In all cases the traditional Breit-Wigner form gets an 
excess in its low-energy wing as is revealed in mass-correlations of the decay 
products. 

Collective flows of final particles present another interesting effect. 
These flows are generated by the internal pressure in the quark-gluon medium. 
Their shapes are determined by spatial characteristics of initial overlap vol- 
umes and by interaction dynamics. The azimuthal angle distribution of emit- 
ted particles can be written as 



where \l/ r denotes the reaction (event) plane determined, independently for 
each component, from the flow itself: v n =< cos[n(0 — ^ r )] > with average 
over all particles in all events as indicated by < >. The first harmonic coef- 
ficient of the Fourier expansion of azimuthal distributions, v±, describes the 
isotropic radial flow, and the second coefficient, i>2, characterizes the elliptic 
flow. To the anisotropic flow there corresponds the azimuthal asymmetry of 
particle distribution with respect to the reaction plane defined by the beam 
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direction and the vector of impact parameter. The higher harmonics are also 
studied. The "non-flow" contributions due to jets, ridge, resonance decays, 
HBT effects are usually subtracted. 

The radial flow manifests itself in flattening (in AA collisions compared 
to the pp ones) of the spectra at small transverse kinetic energy which is 
stronger for heavier particles. It scales with energy and the system size. 
This effect is more pronounced in central collisions between equal spherical 
nuclei. The change of sign of radial flow of protons (the so-called "wiggle") as 
a function of rapidity for the most peripheral events at SPS energies has been 
noticed. No wiggle is observed at RHIC. The coefficient Vi stays negative 
(antiflow) and steadily increasing for all pseudorapidities. This shows the 
early start of collective effects. 

The most intensively studied elliptic flow is related, first of all, to the 
non-centrality of collisions (the spatial eccentricity of the reaction zone) and 
to the rescattering of partons. The elliptic shape of the overlap region gives 
rise to different pressure in different directions and, therefore, to non-zero v 2 . 
Thus, the coefficient v 2 is sensitive to the system evolution at very early times. 
The integrated elliptic flow increases with energy. This increase is stronger 
in more central collisions than in mid-central ones while the absolute values 
are smaller. The value of v 2 depends on pseudorapidity with maximum at 
i] = and increases with energy as well. It drops about a factor of two 
around rj = 3. The value of v 2 is smallest for more central collisions. The 
dependence on the transverse momentum has been measured up to pt = 12 
GeV. It shows the maximum about px = 3 GeV with v 2 becoming as large 
as 0.2 and then some saturation or slow decrease. At p T < 2 GeV v 2 is 
characterized by the scaling with the transverse kinetic energy m T — m and at 
2 < pt < 4 GeV - with the number of constituent quarks. The relation of v 2 
to the spatial eccentricity of an event depends on the number of participants 
(system size). The coefficient v 2 has been measured for different particle 
species. The baryonic elliptic flow saturates at somewhat higher pt and 
larger v 2 values than for mesons. The value of v 2 is close to for J ftp in 
Au-Au collisions at 200 GeV. 

Enhancement of the yield of various particle species at high px in AA 
collisions as compared to the very strong suppression of pions has been ob- 
served. Thus the hadrochemical compositions in AA and pp collisions are 
different. This is considered as a manifestation of the quark-gluon plasma 
(or, more generally, of the prehadronic state) effects. Neutral pions and 
7/-mesons at high px > 2 GeV are strongly (about 5 times) suppressed in 
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central AA collisions relative to pp processes. The suppression of <fi mesons 
is the same at high p t > 5 GeV but smaller in the intermediate 2 < < 5 
GeV range. No suppression is seen for photons in this kinematic range. The 
proton production is enhanced at intermediate pt ~ 2 — 5 GeV so that the 
ratio p/ix becomes close to I. Quite surprising is the fact that in AA at RHIC 
antiprotons are as abundant as negative pions for p T > 2 GeV. That is prob- 
ably related to the collective transverse flow. In particle collisions the strange 
particles production at high pt is strongly suppressed compared to pions (3 
- 4 times) but it is enhanced by about the factor of two in AA processes so 
that the suppression factor becomes about 0.6 - 0.7. The enhancement is 
larger for particles with higher strangeness content such as S, f2. The overall 
fraction of strange particles is about twice as high in heavy-ion collisions as 
compared to the elementary particle collisions. The similar relatively small 
suppression factor was measured for charmed quarks (J/ip)- The upper limit 
for bottom (T) suppression of 0.64 was recently obtained. A strong sup- 
pression of the high-pr electrons originating from the decays of charm and 
beauty hadrons indicates the same value of suppression for heavy and light 
quarks in the quark-gluon medium. This shows that the initial stage of cre- 
ation of highly virtual quarks is more important than their subsequent mass 
difference. The change in the hadrochemical content of collision products 
with high p T may reflect a difference in in-medium interactions of the species 
and asks for taking into account the collective (non-local) properties of the 
medium (with possible diffusion before hadronization). 

Remarkably enough, at very low pr < 0.5 GeV of secondary particles 
their yields seem to be the same in pp, AA and even e + e~ processes (albeit 
different for different species) and independent of the primary energy as is 
expected from universality of coherent processes for long-wavelength gluons. 

The above statements can be illustrated by numerous figures which can 
be found in experimental presentations and more specialized review papers. 
We will only show some of them to elucidate theoretical argumentation. That 
would allow to keep the reasonable size of this article. 

All these observations are closely related to the general physical process 
of the energy loss of partons during nuclei collisions which is to the large 
extent the topic to which the present review is devoted to. 
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3 The microscopic description of the quark- 
gluon medium 

The microscopic QCD description of all processes with thousands of particles 
produced is surely impractical in the whole phase space region. Such a de- 
scription may be applicable to the relatively rare high-p-r subprocesses where 
the QCD coupling constant becomes smaller. For soft partons and hadrons 
the rescattering and hadronization effects are too strong to leave many traces 
of the primary inelastic process. The macroscopic approach (including the 
statistical one) is more suitable for description of purely collective medium 
effects. 

One can nevertheless follow general evolution of fields during heavy ion 
collisions and reveal typical correlation patterns by applying main principles 
of QCD. In this context, the present paradigm, which we describe below, pic- 
tures the transition from the Color Glass Condensate (CGC) to Glasma even- 
tually evolving into Quark Gluon Plasma (QGP) and subsequent hadroniza- 
tion . 

The energy loss of partons in the quark-gluon medium is the main source 
of experimental information about its properties during these stages. As in 
electrodynamics, it may be separated in the two categories. 

The loss due to the change of the velocity vector of a parton such as elas- 
tic scattering, bremsstrahlung and synchrotron radiation is usually treated 
microscopically. All these phenomena result from the short-range response of 
the parton to the impact of the matter fields. Elastic scattering does change 
the parton energy due to the recoil effect and deflects it and thus changes 
the energy flow in the initial direction. At high energies this process is less 
probable than the emission of gluons (bremsstrahlung) due to the nearby 
collisions with the matter constituents. The medium structure imposes some 
restrictions on the coherence properties of the radiation process and on the 
effective radiation length (in analogy to the Ter-Mikhaelyan and Landau- 
Pomeranchuk effects in electrodynamics). In the general case, the intensity 
of the radiation depends on the relation of the path length of the parton in 
the medium, its mean free path (depending on the distance between the scat- 
tering centers and on the cross section) and the formation length of emitted 
radiation. The synchrotron radiation of gluons induced by the curvature of 
parton trajectory in chromomagnetic fields may become important for strong 
enough fields. These processes are considered in the section 3. 
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A principally different source of energy losses is connected with the medium 
polarization by the propagating parton. Here one speaks about the collec- 
tive response of the medium related to the non-perturbative long-range in- 
terconnection of its constituents described by its chromopermittivity. This is 
treated macroscopically within the in-medium QCD. The corresponding ef- 
fects are analogous to the Cherenkov radiation, the wake and the transition 
radiation. In these processes one can, in the first approximation, neglect the 
change of the velocity vector of a parton. The macroscopical aspect of the 
problem is sometimes ignored in the review papers and we will discuss it in 
more details in section 4. 

Let us note that the total energy loss is the sum of the losses due to 
short-range and long-range interactions. 

3.1 Color Glass Condensate, Glasma, Quark-Gluon Plasma 

Is it possible to work out a quantitative description of multiparticle produc- 
tion in hadronic and nuclear collisions within QCD? An ambitious attempt 
of answering this question is provided by the paradigm of Color Glass Con- 
densate (CGC), see e.g. the reviews [H |2j |3j H] . In the words of one of the 
principal creators of this paradigm, "Color Glass Condensate is a universal 
form of matter which controls the high energy limit of all strong interac- 
tion processes and is the part of the hadron wavefunction important at such 
energies" [5]. 

The strength of the CGC approach to the physics of high energy scat- 
tering is in providing a unifying concept for such important but previously 
separately existing notions as soft (wee) and constituent partons, stringy and 
jetty particle production, structure functions and limiting fragmentation, etc. 

Multiparticle production in QCD is generated by low x partons in the 
nuclear wavefunctions (CGC). Chromoelectric and chromomagnetic fields in 
the medium under consideration are the strongest fields in nature. When two 
"pancakes" of CGC collide, they form matter with very high density called 
Glasma. It preceeds the Quark-Gluon Plasma (QGP). Glasma consists of 
color flux tubes localized in the transverse plane and stretching between the 
valence color degrees of freedom. These tubes generate long-range correla- 
tions responsible for flat distributions in pseudorapidity Arj. The existence 
of long-range (pseudo)rapidity correlations is related to the early times of 
Glasma formation in nuclear collisions r < r out exp(— \y a — j/&|/2), where r out 
is the proper time of production of the particles a and b with the rapidities 
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y a and yj,. The radial flow collimates particles in the direction of the flow 
and generates azimuthal correlations but does not affect the distribution over 
A77. In particular, the ridge could be explained as the combination of the 
initial state long-range rapidity correlations and the final state radial flow. 
Therefore combined elements of the both micro-and macro- approaches are 
necessary to explain the ridge phenomenon. 

Color Glass Condensate is an effective theory describing hadronic and nu- 
clear wavefunctions in the regime where the main contribution to the cross- 
section is coming from their multiparton Fock components, i.e. in essence 
from the dense gluon matter. The key idea behind the physical picture of 
CGC can be formulated as follows. Let us consider a fast hadron (nucleus) 
moving along the z axis with the big longitudinal light-cone momentum P + . 
The modes in the projectile are classified by their relative yields of longitudi- 
nal momenta x = p + /P + . Then, from uncertainty relation, the characteristic 
longitudinal extension of these modes is Ax~ ~ 1/ (xP + ) and their character- 
istic lifetime is Ax + ~ — 2xP + /m 2 ± Q where m\ is some characteristic 
scale of transverse momentum. We see that the hard constituent modes, for 
which x ~ 1, are localized at small distances and have large lifetimes, while 
the soft wee modes, for which x <C 1, are spread over large longitudinal 
distances and are short-lived. From the data on structure functions from 
deep inelastic scattering we know that the number of wee partons, most im- 
portantly gluons, is rapidly growing with 1/x, xG(x,Q 2 ) ~ (l/x) K with 
K ~ 0.3, so that these modes are characterized by large occupation numbers. 
For such modes it is natural to replace the language of particles by that 
of fields, and, in turn, for the coherent multiparticle states the appropriate 
description is that in terms of classical fields. 

These considerations lie at the origin of the McLerran-Venugopalan model 
of nuclear light-cone wavefunctions P, [7J E] in which, in the framework of 
QCD, the constituent and wee partons are integrated in the unifying picture 
in which strongly localized slowly evolving constituent hard modes are con- 
sidered as the source for the delocalized soft wee gluon field for which one has 
to take into account nonlinear effects leading, in particular, to the saturation 
of the gluon distribution, see below. 

Let us introduce the variables to be used in this section. Most of the 
formulae are written in terms of light cone coordinates. For the four-vector 

4 The light-cone coordinates x + and x~ are defined below in equation fcih 
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x M = (t,x±,z) the definition of the ±-components reads 

The above-described momentum components are described analogously. We 
shall also make an extensive use of the (r, 77) coordinate system 



(4) 



r = V 2x + x~, r] = - In f — j , 
2 \x-J 

where r is the proper time and r\ is the spatial rapidity 

3.1.1 McLerran-Venugopalan model. Saturation. 

To illustrate the machinery of the McLerran-Venugopalan model, let us con- 
sider a calculation of the momentum space density of gluon modes 

dN -> -, 2k + 

— = (a?(x + ,k)al(x + ,k)) = j^(A\{k,x + )A\{-k,x + )), (5) 
from which one can calculate the observable gluon structure function 



xG(x, Q 2 ) = J d 2 k ± k 



dN 



dk+d 2 kt 



(6) 

x=k+/P+ 



As we are interested in the high density regime, the gauge potentials A l a (k, x + ) 
in (j5]) correspond to soft wee gluon fields. 

To calculate the momentum space gluon density in ^ one needs to spec- 
ify the gluon fields A l a (k,x + ) and the nature of averaging. The first task is 
achieved by assuming that A l a (k,x + ) satisfy the classical Yang-Mills equa- 
tions with the source originating from hard constituent modes 

[D ll ,F^=6" + p 1 (pc ± ,x-), (7) 

where F^ u = [D^, D v ] is the non-abelian field strength and _D M = — igA^. 
The solution of ^ is conveniently written in the light-cone gauge A + = 0. It 
is easy to verify that for the static x + -independent solution one has A~ = 
so that the solution of ^ is purely transverse: 

A*(x ± ,a;-) = -C/(x ± ,a;-)9 i [/ t (x ± ,x-) ) (8) 
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where 

U(x±, x~) = P exp i^ig J dy~a(~x±, x~) \ , (9) 

and a(x_L, x~) = — p(xj_, x~)/V\ is the solution of Q in the covariant gauge. 
The solution g,(|9|) is a nonlinear functional of the color charge density 
p(x±,x~) and, being exact, is nevertheless formal. 

The second crucial step is related to the above-mentioned specification of 
the nature of averaging in ^ . The idea of [SI [7J [8] was that for a large nucleus 
and the transverse resolution of the probe | Axj_| <C 1/Aq CD , where Aqcd 
is the fundamental scale of strong interactions, Aq CD ~ 100 MeV, the color 
density p(xj_,x~) is, on the event-by-event basis, essentially a random func- 
tion. In the simplest setting where the hard modes are constituent quarks 
the static color density is just the total color charge of all constituent quarks, 
that happen to lie in the tube under consideration. An averaging procedure 
in ([5]) is thus specified by some functional W\+ [p], where A + is the scale sep- 
arating hard modes from the soft ones. A simple but quite realistic averaging 
is provided by the Gaussian ensemble such that 

(p a (x ± ,x-)p\y ± ,y-)} WA+ = g 2 ^ A 5 ab 6 2 (x ± -y ± )5(x- -y-)). 

In the simplest case where the sources are constituent quarks the factor p\ 
setting the scale for the transverse density of color charge is simply p 2 A = 
g 2 A/{2nR 2 A ) and is thus proportional to A 1 ^ (see e.g. a detailed derivation in 
[T]). The color charge density pa is directly related to the physical saturation 
scale Q s to be defined below, Q s « 0.6 g 2 pA [ID] - 

It turns out convenient to present the expression for the gluon density ^ 
resulting after averaging over the stochastic sources p with the above-defined 
Gaussain weight for the gluon density in transverse space (p(x, J] 

2 4vr 3 1 d 3 N 
<P{*,*i) ~ N 2- l7T R A d\n(lMd 2 kx (10) 

The characteristic properties of the resulting distribution can conveniently 
be presented by writing down its asymptotic expression^] at k. 2 ± — > and 

5 For a detailed analysis supporting this assumption see [5]- 

6 The convenience of introducing tp(x, Kj_) lies in particular in the fact that the quantum 
evolution of gluon distributions considered below is most naturally written precisely in 
terms of this so-called unintegrated structure function. 

7 For details of this calculation see, e.g., [T]. 
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k]_ -> oo: 

i/2 



(11) 



The most important ingredient of (11 ) is the so-called saturation scale Q s ~ 



A 1 / 3 . The A-dependence of is directly related to that of 

The remarkable phenomenon of flattening of the gluon distribution at 



small transverse momenta seen in (11) is called saturation. We see that sat 



uration is characteristic of the regime in which the gluon density (p ~ l/a s 
is parametrically large. Another crucial observation which is of great signif- 
icance for developing the quantum theory of gluon distributions is that at 
large A the saturation scale becomes large, Q 2 S ^> Aq CD . Since the saturation 
momentum Q s sets the scale for all the transverse momenta in the problem, 
for asymptotically large nucleus the theory is in the weak coupling regime. 

3.1.2 Glasma 

In the paragraph 3.1.1 we have described, in the framework of color glass 
condensate theory, the distribution of gluons inside the nucleus. With this 
knowledge at hand we can consider a problem of particle production in high 
energy nuclear scattering. The main idea of CGC-motivated approach to 
this problem is that, if gluon distributions in nuclei before the collision are 
most naturally described in terms of fields, it makes full sense to describe the 
energy and particle flux created at the early stage of these collisions also in 
terms of gluon field. In the tree-level approximation this amounts to solving 
the Yang-Mills equations with two external currents pi,2(xj_, x~) representing 
color distributions of the two incident nuclei propagating along x = 0: 

[D„F^] = 6» + Pl (x ± ,x-) + ^ + p 2 (x ± ,x-). (12) 

If we are interested in calculation of the gluon production in the central re- 
gion, where one can assume that their rapidity distribution is approximately 



flat, it is natural to look for the boost-invariant solution of (12) in the form 
[II]: 

A* = 6{-x + )6{x-)A\ 1) + 6(x + )6(-x-)A\ 2) + 6(x + )6(x-)A\ 3) 
A v = 0(x + )0(x-)^ 3) , 
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Here A\ 2 ^ are the gluon fields of the two incident nuclei before the collision 
described by (8) and A\^ and Aj 3 \ are the gluon fields created after the 
collision that is assumed to take place at r = 0. Taken at the proper time 
t = + , these fields describe the initial conditions for subsequent evolution 
of dense matter created in high energy nuclear collisions. 

It is straightforward to verify that to ensure the matching of the solution 



of (12) at r = the following relations should take place: 



^(3)|t=0 — -A(l) +^4( 2 )) 



4)' r=0 - 17 [4)'4)] • ( 13 ) 



A crucial observation [T2| IT3] following from the boundary conditions (13) is 
that at r = + there exist only longitudinal fields 

E z = iglA^Afa], 



B z = ige 



Kl 



4> 4) > ( 14 ) 



in which the energy created in the collision is stored. Thus there takes place 
a full transformation of the structure of the solution, i.e. the transition 
from the two transverse fields A\ x 2 y described in the framework of CGC at 
t = 0~, to the purely longitudinal chromoelecric and chromomagnetic fields 
immediately after the collision at r = + called glasma. Let us stress that, 



as obvious from (14), the effect is purely non-abelian. As the only scale in 
the problem is the saturation scale Q s , the correlation length in the impact 
parameter plane is R± ~ 1/Q S - Thus, the physics is that of the longitudinal 
fields characterized by the transverse correlation radius R±, i.e. with a set 
of the chromoelectric and chromomagnetic cylindrical tubes. It is at this 
point where the CGC description matches the old paradigm of high energy 
hadronic collisions, the Lund string model. An important new ingredient, 
however, is the presence of the chromomagnetic flux tubes. The very fact 
that E ■ B ^ opens up a principal possibility of the local CP violation [13] . 

In the course of subsequent evolution the initial configuration of longitu- 
dinal fields is undergoing a partial transformation into that also containing 
transverse field and at r ~ 1/Q S the energy stored in longitudinal and trans- 
verse modes is approximately equal [13]. 

The expressions for the number of initially produced gluons and their 
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transverse energy can be written from dimensional reasons^] 



dN 
dr] 



c N 



ri=0 



*R 2 A Q\ 

a s (Qs) 



dE\ 



dr] 



r)=0 



a s (Q s ) 



(15) 



where the coefficients ~ 0.1/47T and ce ~ 0.05/47T, (see e.g. [15] and 



references therein), are determined from numerical solutions of (12). 

In particular, using the value of Q s determined from the analysis of deep 
inelastic scattering at HERA for x e g ~ 0.01 as appropriate for nuclear scat- 
tering at RHIC, Q S (RHIC) ~ 1.2 GeV, one gets from (ph the rapidity 



density of the number of gluons and of the transverse energy: 



dN 
dr] 



1100 



dE\ 



77=0 



dr] 



500 GeV, 



(16) 



r)=0 



which agrees with the experimental data. This can be interpreted as an in- 
dication of almost perfectly ideal isentropic evolution from the initial glasma 
phase, through the intermediate QGP stage, to that of the flux of final 
hadrons. 



3.1.3 Glasma instabilities 

In the section 3.1.2 we have considered a boost-invariant solution of the 



classical Yang-Mills equations (12) describing the conversion of two initial 
gluon fields inside the incident nuclei into the physical fields created in the 
collision. The latter can be described as a set of chromoelecetric and chromo- 
magnetic flux tubes having the radii of order of inverse saturation momentum 
Q s stretched between the receding sources. 

In general, one expects that classical solutions of Yang-Mills equations are 
unstable with respect to quantum fluctuations. For example, constant chro- 
moelectric fields are unstable with respect to Schwinger gluon pair production 
and static chromomagnetic fields are characterized by the Nielsen- Olesen in- 
stability. A tree-level iterative calculation of gluon production in collisions of 
the non-abelian charges in [TO] has also demonstrated the instability of the 
approximate tree-level solution with respect to higher order corrections. 

Thus it is probably not too surprising that the evolving glasma fields are 
unstable with respect to rapidity-dependent fluctuations. Such fluctuations 



8 Because of the ultraviolet divergence at r = these quantities can meaningfully be 
computed only at nonzero r > [14] . 
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are inevitable when one takes into account quantum fluctuations. An explicit 
calculation of the wavefunction of two colliding nuclei demonstrating this 
effect was done in [TJj. 

In describing the corresponding instabilities it is convenient to consider 
two regimes: 

• Small proper times r < 1/Q S . This is a regime in which instabilities 
are those related to fields. 

• Intermediate times r > 1/Q S . This regime can be described in dual 
terms of fields and particle kinetics, see [T8] . 

The analysis of these instabilities was brought to a new level in the all- 
order numerical calculations made in [T$t |2"0"] in which it was shown that the 
amplitudes of the soft modes grow like exp(y / r). An analytical analysis of 
instabilities with respect to quantum fluctuations in the flux tube picture 
was performed in [21]. 

The physical origin of these instabilities is still under debate. 

The closest known analogy is the non-abelian Weibel instability [221 [231 
[21] with plasma instabilities of both electric and magnetic type developing 
in non-equilibrium plasmas with inhomogeneous distributions of particles 
in momentum space. The analogy to glasma is very natural because, as 
described in the previous paragraph, the initially created configuration of 
glasma fields is maximally anisotropic (the initial transverse momentum is 
equal to zero). 

The analysis carried out in [21] have lead to the complete revision of 
the previously existing picture of particle production in nuclear collisions 
given in the framework of the so-called "bottom-up" scenario of thermal- 
ization of the glue initially freed at the momentum scale Q s described in 
|25j . The reason for revision was that instabilities are much more effec- 
tive in producing soft gluons and driving the system towards isotropy faster 
than the bremsstrahlung mechanism considered in [25]. The new versions 
of the "bottom-up" scenario include, in particular, mechanisms of cascade 
turbulent-like multiplication of gluons [2SJ [23 I2H] • 

3.1.4 Quantum evolution. High-energy factorization. 

The McLerran-Venugopalan model is by construction adequate for moderate 
values of fractional light-cone momentum x = k + /P + . The quantum cor- 
rections to the tree-level correlators like rt5J) involve soft gluon emission and 
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virtual gluon loops. Both effects bring in contributions which paramerically 
are of the order of 0(a s log(l/x). As physical quantities can not depend 
on the arbitrary separation scale of the soft and hard modes log(A + /P + ), 
there arises a Wilsonian renormalization group resumming such contribu- 
tions. For dense gluon systems this resummation of small-x enhanced terms 
has to include contributions to all orders in the gluon density. The loga- 
rithmic nature of quantum corrections makes it natural to use the rapidity 
variable Y = log(l/x). In these terms one can, for example, interpret the 
gluon two-point function ^ as an average of the product of classical fields 
taken at some specified rapidity scale: 

< AA > Y = J VpW Y [p] A cl (p) A d .(p), (17) 

The Wilsonian RG equation (the Jalilian-Marian-Iancu-McLerran-Leonidov- 
Kovner (JIMWLK) evolution equation (29J, [301 |3T] ) is most naturally written 
as an equation for the functional Wy[p]. As shown in [32] it can be cast in 
an elegant Hamiltonian form: 

dW Y [ P } 

dY 

where H[p] is the JIMWLK Hamiltonian 

H[P] = o / X(x-L,y±) T- f u 7 x + / <t(xl) j-j-^-. (19) 
2./x x ,y x 6p{x. ± )5p(y±) J x± op(x±) 



H[p] W Y [p], (18) 



In equation (19) the functions x( x ± ; y±) an d °"( x ±) are explicitly calculable 
[29l I30| l3Tj and correspond to the effects due to emission of real and virtual 
gluons respectively. 

A remarkable identity relating the real and virtual kernels x(xj_, y_i_) and 



<t(x_i_) discovered in [32] allows to rewrite (19) in terms of the real kernel 
X(x±,y±) only: 

"LPS & 
U[p] = - / -— - x(x_L,y_L) — - (20) 

The explicit expression for x( x ±;y±) involves the path-ordered exponentials 
U(x±) defined in equation ([9]): 

[ d 2 z ± (x ± - z ± )(y_L - z_ 

x( x ±,y±) 



4tt 3 (x ± -z ± ) 2 (y ± -z ± ) 2 
[(1 - f/t (x±)[/(z± )) (i _ rf( x± )U(y ± ))] (21) 
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Until recently the usage of the JIMWLK evolution equations and their 
simplified versions, the BK equations [331 EI] corresponding to retaining only 



the two-point functions in the infinite hierarchy generated by (18), was re- 
stricted to one-source problems like deep inelastic scattering at high energies 
and quantities like gluon structure functions. In the major recent develop- 
ment [351 ESI EZ] the formalism allowing systematic computation of quantum 
corrections to the two-source glasma problem for inclusive gluon production 
was developed. In particular, a generalized factorization property of the 
two-source problem in high energy QCD was derived in which both incoming 



fluxes are described by the solutions of the JIMWLK equation ( 18 ) charac- 
terized by the corresponding saturation scales. This approach allows to cal- 
culate the inclusive gluon production cross section. Schematically the leading 
logarithmic correction to such observable O can be written as [351 E3 137] : 

AO = [log(l/x!)Hi + \og(l/x 2 )H 2 } O (22) 

where \og(l/ x\$) are the rapidity shifts taken into account in computing the 
quantum corrections and Hip are the JIMWLK Hamiltonians of the incident 
nuclei. 



3.2 Jet quenching and parton energy loss in dense non- 
abelian medium 

As has been already mentioned in the section 2, one of the most impres- 
sive results of the experiments at RHIC is a dramatic decrease in the yield 
of particles with large transverse momentum - a phenomenon directly re- 
lated to the so-called jet quenching in dense QCD matter. The underlying 
phenomenon is, of course, more general than just the above-mentioned at- 
tenuation of inclusive spectra of transverse momentum but rather refers to 
the specific medium-induced changes in the physics of multiplarticle produc- 
tion. The main focus of jet quenching studies is on looking for modifications 
of standard patterns of processes at large transverse momenta (such as jet 
production in nuclear collisions) as compared to the results for pp collisions. 
The jetty structure of final state in events characterized by large transverse 
energies is one of the generic features of particle production in high energy 
hadron collisions. The properties of these jets, i.e. the collimated fluxes of 
hadrons, are very well described, especially in the e + e~-annihilation, by per- 
turbative QCD under assumption of the local parton-hadron duality. This 
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refers both to the cross-section and topology of jetty events and the proper- 
ties of the individual jets such as intrajet particle distributions, multiplicity 
distributions, etc. Thus QCD jet physics is a very natural place to look 
for medium-induced modifications in nuclear collisions - in the hope that 
one can achieve reasonable theoretical and experimental control on its rel- 
evant characteristics. For the recent reviews on theoretical approaches to 
describing jet quenching and their comparison with experimental data see 

e.g. [3H1 EHl SQl E21 SSI SS]. 

Schematically one can break the process of jet formation into three stages: 

• Production of an initial high energy parton. At asymptotically high 
energies a relevant theoretical description is given by the standard for- 
malism of collinear factorization 

• Creation of final state high energy particle is generically accompanied 
by its radiation. Technically it is convenient to ascribe to the initially 
produced parton some virtuality Ql which is shaken off by subsequent 
perturbative radiation of gluons and quark-antiquark pairs until reach- 
ing some cutoff invariant mass Q\ at which the perturbative description 
is no longer applicable. 

• Hadronization of partonic configuration created in the course of QCD 
jet evolution 

Let us consider, for example, the simplest inclusive distribution (fragmen- 
tation function) D h {x = Eh/Eo, Qq) describing a probability for a hadron in 
the jet originated by the partor^] with the energy E and virtuality Q to 
have the energy Eh = xEq. At sufficiently large one can write 

D h (x, Ql) = I dz V(z, Ql Q 2 Q )D h (x, Q\) (23) 

J X 

where V a b{z, Q\, Ql) is the perturbatively calculable probability of produc- 
ing a parton with the fractional energy z = E/Eq and virtuality Q\ in 
the jet originated by the parton with the energy E and virtuality Qq and 
D {x = Eh/E,Ql) is the nonperturbative fragmentation function describ- 
ing the corresponding parton-hadron conversion. Probabilistic description of 

9 In this section we will omit all parton subscripts and, for simplicity, consider the purely 
gluonic cascade. 
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parton-hadron transformation is possible due to collinear factorization prop- 
erty of QCD that ensures, at large enough Qq, the absence of interference 
between perturbative evolution of the initial parton and the final nonper- 
turbative process of color blanching and formation of the final hadron. The 
probability V(z, Q\, Qq) refers to a QCD cascade converting the initial par- 
ton into a spray of partons with different fractional energies {z} and the same 
final virtuality Q\. In the in-vacuum QCD the evolution of the fragmentation 
function D^(x, Qq) with Qq is described by the standard Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) evolution equations [4*6l Wf\ |4"8]. 

At large Qq the formation of the initial parton takes place at very small 
time and, in the first approximation, the medium effects on this initial hard 
process can be neglected. However, in the in-medium QCD both the proba- 
bilistic kernel V(z, Q\, Qq) and the nonperturbative fragmentation function 
D h (x = Eh/E,Ql) depend, generally speaking, on the properties of the 
medium. The dominant opinion in the current literature is that all the rel- 
evant medium-induced effects take place at the perturbative stage and refer 
to V(z,Ql,Ql) leaving D h (x = E h /E, Q\) untouched Following this as- 
sumption, one can introduce the medium-modified fragmentation function 
D h a {med \x,Ql) 

D h J™ d \x,Ql)= I dz(V ab (z,Ql,Ql)+AV ab (z,QlQl\ {//})) D%(x,Q 2 h ), 

J X 

(24) 

where AV a b(z,QhiQo\ {/*}) * s the medium-induced change of the in-vacuum 
probability that depends on the properties of the medium under consideration 
that are characterized by the set of parameters {/i}. 

To calculate AV(z, Q 2 h , Ql\ {/i}) one has to understand the mechanisms of 
medium-induced effects such as medium-induced radiation and elastic energy 
loss experienced by a parton propagating in the dense non-abelian medium. 
This is discussed in the paragraph 3.2.1. 

Let us note that both of the two above-described medium-induced mech- 
anisms are well studied in the abelian case in the framework of the physics 
of the electromagnetic showers in matter (see e.g. [12]). An important ingre- 
dient of the physics of electromagnetic showers is accounting for ionization 
losses leading to specific distortions of the pattern resulting from medium- 
induced bremsstrahlung and pair creation [39]. The modifications of QCD 

10 Experimental data which favor it were mentioned in section 2. 
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cascades taking into account " pionization" , i.e. the interaction with addi- 
tional strongly interacting degrees of freedom, were analytically studied in 
[501 EE E21 E3] for the case of non-angular ordered cascades and in [5U [55] 
for the angular-ordered ones. 



3.2.1 Parton energy loss 

In the previous paragraph we have mentioned that the process of transfor- 
mation of the initial highly energetic highly virtual parton into the spray 
of final hadrons is, in some approximation, that of a perturbative cascad- 
ing emission of new quanta. Thus we are dealing, e.g., with a problem of 
multiple gluon emission. To compute generic in-medium corrections to the 
in-vacuum picture one has, in principle, to be able to calculate cross sections 
for multiple medium-induced gluon emission. As of now this problem is not 
solved. At the current level of the theory of the medium-induced radiation 
one can compute the corresponding cross sections for one-gluon emission in 
the first order in the strong coupling constant but, under certain simplifying 
assumptions, to all orders in the interaction with the medium. The resulting 
spectrum of one-gluon emission can be written as a sum of the vacuum and 
medium-induced contributions 

d jtot rfJ vac Jjmed 



where 



dwd 2 k± dwd 2 \s.± dwd 2 \t 

dP ac a, N, 



+ (25) 



dwd 2 k i 2it 



K™ c (z), (26) 



and K vac (z) = — -—2+z(l — z) is the standard DGLAP gluonic splitting 

z{l - z) 

function. The general case with multiple gluon emission is then considered 
through constructing a probabilistic branching process obtained by iterating 



the one-gluon emission described by the first-order expression (25). The 
corresponding theoretical schemes used for constructing such cascades will 
be discussed in the paragraph 3.2.2. Below we shall discuss the methods of 



calculating the second term in the right-hand side of (25) to the accuracy 
0(a s ). 

The physical problem at hand is the computation of the leading contri- 
bution to the medium-induced energy loss AE of the parton produced inside 
the medium on its way from the hot dense fireball. Generically the energy 
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loss AE depends on the medium thickness L and the bulk characteristics of 
the medium and is described by some probability distribution V(AE). 
More precisely, the losses are determined by the following factors: 

• Probability of an elementary event leading to energy loss characterized 
by the opacity N = L/X, where A is the mean free path of the parton 
in the medium under consideration. For the particle with integrated 
particle-medium interaction cross section a and the medium density p 
the mean free path can be estimated as A ~ l/(pa). 

• The intensity of the impact or scattering power of the medium charac- 
terized by the transport coefficient q = {p±)/X, where {p 2 ± ) is the aver- 
age transverse momentum squared that the propagating particle gets 
from the elementary act of collision. In the thermal medium q = m 2 D /X 
where mo is the Debye mass. 

There exists a general consensus, supported by extensive model calcula- 
tions, that the radiative medium-induced energy loss is the dominant one ( see 
e.g. [56]). Theoretical work on studying the properties of medium-induced 
gluon radiation was carried out by several groups [571 E3 ES EH ED |62j E21 



EligiiaEBEElIisiTjB^ The 



approaches of the above-listed groups differ in details of treating the relevant 
kinematics of the radiation and the details of describing the medium under 
consideration. In the present review we will sketch, following jlOl H3J EH] , the 
main steps of the approach developed in [671 ESI EH EQl CD G21 E3J [7J1 [75] . The 
two key ideas allowing to carry out the calculation of the medium-induced 
radiation in the form discussed below are the separation of transverse and 
longitudinal degrees of freedom [SJ] and the usage of the light-cone path 
integral formalism [72J. 

Let us start with the fundamental idea that, at variance with QED, in 
QCD one can meaningfully consider radiation of gluons in the eikonal approx- 
imation. This radiation takes place due to decoherence of the initial coherent 
parton flux by soft color rotations caused by the fields of the medium (see 
a clear discussion in [301 SI])- The basic objects of the eikonal approxima- 
tion to QCD are the Wilson lines. Staying in the A + = gauge used in 
the previous section, we get the following expression for the Wilson line in 
question: 




(27) 



26 



where x_|_ is the coordinate of the incident parton in the impact parameter 
plane which in the eikonal approximation stays constant and A~(x + , xj_) are 



the target fields. The Wilson line (27) fully determines the S-matrix element 
of the eikonal scattering 

S{p ^p')~5 (p' + -p + ) J dx ± e- ix ^ p t- p ^W(x ± ). (28) 

The cross section is therefore fully specified by the average of the product of 
Wilson lines taken at different points in the impact parameter plane 

\S(p^p')\ 2 ^jW*(x.±)W(y ± )). (29) 

To elucidate the physical meaning of the average of the product of Wilson 



lines in (29) let us consider the simplest and at the same time most popular 
model of the medium in which the medium is described as a collection of 
static sources of field a~(x. ±i ) located at points {xf} and thus characterized 



by the density n(x + ) = ^2S(x + — xf). For such a medium one gets for (29 ): 
i r (Wrf(x 1 )W(y_ L )) ^expj-^ J da; + n(*Xx ± -y ± )}, (30) 
where 

^(xi -y±) = J |^Tr| a-(qj| 2 (l - e ^~^) (31) 

is the cross section of scattering of color dipole on the target under con- 
sideration. Let us stress that the dipole cross section appears only in the 
matrix element squared (or, equivalently, the cut diagram); the problem 
at hand is that of scattering a single color probe. In the soft scatter- 
ing approximation the dipole scattering cross section is simply expressed 
as <t(x_i_ — y_i_) ^ C(x_i_ — y_i_) 2 , so that 



L I (M/t( x± )W(y ± )) ~ exp j^-L J dx + q(x + )(x ± - y ± ) 2 | 



(32) 



where g(0 = 2^2Cn(C). Here q is nothing else but the above-introduced 
scattering power of the medium. 

If one is interested in going beyond the soft scattering approximations, it 
is necessary to take into account the momentum transfer from the medium. 
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The corresponding elegant generalization of ( 29 ) reads 



S(p ->p') ~ 5 (p '+ - p + ) p + J dx ± e^ Xx)(p ^ p±) 

dr±(x~ 



Vr±(x ) exp < i 



dx 



dx^ 



W(r ± ). (33) 



In (33) the momentum transfer from the medium is described by the motion 



of an effective nonrelativistic particle (a standard situation in the light-cone 
quantization) in the impact parameter plane. 



The formula (33) allows to calculate the spectrum of medium-induced 



gluon radiation in the medium of longitudinal extent L + that, in particular, 
takes into account the interference of in- vacuum and in-medium contributions 
[SZj (see also [7S]): 



dl 



dk+d 2 k 



a s C R 1 a s C R 

7T 2 k 2 



2n) 2 k+ 




2Re^ + dxf j d 2 x ± e~ ik±xx 







dx+e rU*+ «»«M*J-) 
k + .L+ ' dy± (9xj_ 



/C(xj_, x\ | 0, x\ 







7C(x ± ,L+| 0,x+ 



k 2 ± dy ± ' 

Here /C(rj_(x^")| rj_(a;^)) is the Green function 



(34) 



JC(r±(x^)\ r ± (xt)) = / Pr ± exp 



(35) 



In the soft scattering approximation the Green function (35) simplifies to 

■ P 



]C(r ± (x2) \ r ± (xf)) = / Vr ± exp < i 



di [ii + i- 



• (36) 



2 ./, V~~ L 2y/2p+ 
The expansion in opacity or, equivalently, in the density of the medium 



is obtained by expanding (34) in the cross-section a. 



The physics behind the full expression ( 34 ) can be illustrated by consid- 

- 2u/k\. 



ering a soft gluon emission. The formation time of a gluon is Tf 
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This time should be sufficiently long to ensure the phase change A$ ~ 1 
sufficient for decoherence and thus radiation. In more details, 



A$ = (— ) 

Tform 



k\ 

2w* 



■Az 



qL 
2u 



L 



OJ 



where 



(37) 



(38) 



In the limit oj <C oj c the medium acts coherently reducing the intensity of 
radiation while that of oj 3> oj c corresponds to incoherent scattering. 



OJ 



OJ 



dl 

du 
dl 

du 



a. 



qL 2 



OJ 



O'.s- 



qL 2 



OJ 



In the soft limit u<u c one has for the average energy loss 

dl 



(AE) 



dujuj- 



doo 



qlS. 



(39) 



(40) 



For some time there existed an opinion that the quadratic dependence on the 
medium thickness L is a specifically non-Abelian effect. Recently, however, 
in [25] (see also [86J ) , it was shown that the effect is generic and is valid both 
in QED and QCD for the radiation of particle produced inside the medium 
at the initial piece of its trajectory. 



3.2.2 Multiple gluon emission. In-medium QCD cascades 

At the beginning of this paragraph let us consider the simplest version of 
taking into account the multiple medium-induced gluon emission in which 
one neglects the evolution of the timelike QCD shower in virtuality. Here, the 
possibility of multiple emission is taken into account by constructing a simple 
Poissonian picture in which the probability distribution of the cumulative 
energy loss V(AE) reads 



V(AE) = PoJ2~il II / duj ^ * 

n=0 J i=l J 



dl 



mcd 



Po = exp 



— / dudk i 



' duidk± i 

d jmed 1 

dojdk 



j=i 



(41) 
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The Poissonian cascade described by (41 ) can be derived from the full DGLAP 



evolution (see the Appendix in [87]). Knowing V(AE) one can compute the 
distribution V(e) of the fractional energy loss e = AE/Eq of a parton with 
the initial energy E and obtain the equation for the medium-modified frag- 
mentation function originally suggested in 



^ = Jtt- £ v ^ v ^ (i^' Q 2 ) ■ ( 42 ) 

Turning now to the discussion of the complete picture let us stress once 
again that the energy flow in-vacuum high energy processes involving quarks 
and gluons that are characterized by large momentum transfer is dominated 
by well-collimated jets originating from multiple cascading radiation of glu- 
ons and quark- ant iquark pairs. The physical origin of the cascading process 
is the intense gluon radiation of a high energy parton created in a hard sub- 
process. From the operational point of view one can consider this radiation 
as a cascade degradation of the (large) virtuality of the initial parton un- 
til reaching some non-perturbative scale Qh at which preconfined colorless 
clusters and/or QCD strings are formed. 

To make detailed studies of in-vacuum QCD cascades one turns to Monte- 
Carlo simulations. The two most popular versions of these MC generators are 
PYTHIA [EH] and HERWIG (Hadron Emission Reactions With Interfering 
Gluons) [HIS EJ E2] • Both of them implement the key feature of in-vacuum 
QCD cascades, the angular ordering of successive decays that has its origin 
in fine-tuned quantum coherence effects. In the older version of PYTHIA 
the angular-ordering restrictions were implemented by hand while in its later 
versions and in HERWIG they are incorporated by construction by choos- 
ing an appropriate evolution variable. In the presence of the medium there 
appear new mechanisms that change the structure of in-vacuum QCD jets. 
In the in-vacuum cascade the only source of evolution is the initial virtual- 
ity of a jet. In the in-medium case the medium acts as a source of "extra" 
virtuality to the propagating partons resulting in induced radiation and a 
sink of energy due to the energy loss resulting from elastic and inelastic colli- 
sions of propagating partons with the particles in the medium. The first MC 
code adding medium-induced radiation was PYQUEN (Pythia QUENched) 
[93] . Recently there appeared MC codes describing the medium-induced ef- 
fects in QCD cascades such as JEWEL [94J, Q-PYTHIA [95] (both based 
on the mass-ordered version of PYTHIA) and "Q-HERWIG" [HE] based on 
HERWIG. 
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Let us reiterate that the complete picture combining the in-vacuum and 
in-medium evolution is still not available. The main difficulty lies in the 
fact that the medium created in high energy nuclear collisions is necessarily 
of some finite extent, so that the modifications induced by the medium are 
dependent on its characteristics in space-time: one has to place the cascade 
inside the medium, see Fig. 1 in which the origin of the cascade is placed 
at some distance L from the border separating the dense medium from the 
normal vacuum. 




Figure 1: In-medium QCD cascade. 
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Obviously, the larger the distance of parton propagation in the medium, 
the bigger are the medium-induced losses. The key question is therefore to 
choose, for describing the in-medium QCD cascade, the language allowing a 
credible space-time interpretation. As emphasized in [94] . the only scheme 
allowing such interpretation is the old-style PYTHIA virtuality-ordered cas- 
cade where connection to the spatiotemporal pattern is achieved by calcu- 
lating the lifetime r of a virtual parton which, for a parton with the energy 
E and virtuality Q 2 that has been created in the decay of its parent parton 
with the virtuality <3p ar , reads 



r = E[±- 7 ±-). (13) 



Q 2 Ql 



par 



The lifetime of the initial parton is taken to be r- m = E- m /Qf n . For example, 
the time at which the parton with the momentum q\ in Fig. 1 produced 
the final prehadron is 



^ 1 „ / 1 1 \ E Q 
r h = E — + Ei ] 



Q 2 1 \Qi Q 2 J Q 2 



2 



where x = E\/Eq. We have already mentioned that in the absolute majority 
of studies devoted to energy loss in high energy nuclear collisions it is as- 
sumed that hadronization takes place outside of the hot and dense fireball. 



Setting out the clock with the help of the formula (43) helps to address this 
issue quantitatively [97]. For example, for the parton q\ in Fig. 1 it hap- 
pened that Th < L so that both its formation at time Eq/Qq and its decay 
into final prehadron took place inside the medium. To get a feeling for the 
corresponding physics, let us consider, for example, a gluon with the energy 
Eq = 100 GeV equipped with a typical initial invariant mass Qo = 10 GeV. 
This gluon will on the average decay at r ~ E /Ql = 1 GeV" 1 ^ 0.2 fm, 
so that if the vertex in which the initial hard gluon was generated was suffi- 
ciently far (several fermies) from the surface of hot fireball, a significant part 
of cascade vertices will in fact be generated inside the medium. 

In what follows we shall concentrate on the description of the old-style 
PYTHIA virtuality-ordered cascade and its in-medium modifications. The 
cascade evolution goes through a sequence of decays q =>- q± + where 
q = (E, q) is a four-momentum of the parent gluon and qu2) are those of the 
daughter ones. We study the timelike evolution in which cascading reshuffles 
the initial jet virtuality in such a way that at each decay one has q 2 > q\ + q\- 
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The cascading process at a given branch stops when a virtuality of the last 

11 



parton reaches a threshold level Q\, i.e. no decay into two partons with 



minimal virtuality can occur. The key property of the in- vacuum QCD 
cascades, effective angular ordering of subsequent decays due to color coher- 
ence (i.e. ensuring that, e.g., 9 2 < 9 1 in Fig. 1), is not automatically taken 
into account in the mass-ordered evolution and has to be enforced explicitly 
by accepting only those generated decays for which the condition of angular 
ordering does hold. The basic ingredient of the probabilistic formalism of 
QCD cascade evolution is the conditional probability for a parton produced 
at the scale Qf n with the energy E to undergo a subsequent decay at the 
scale Q 2 < Q 2 n into two daughter partons with the energies zE and (1 — z)E 
correspondingly that can be calculated from the Sudakov formfactor 



S(Q 2 \E,QL;Q 2 h )= exp 



Ql rz+(E, kjj Q\) jjtot 



.(E,kl\Ql) dzd 2 k± 



(44) 

where z±(E, k\\ Ql) define the kinematic limits available for the parton 
decay 

i / n nn / n 2 \\ 

(45) 




and where the full differential probability of gluon emission dI tot j dzd 2 kj_ was 



defined in (25). 



When considering the properties of the in-medium QCD cascades one has 
to deal with several mechanisms modifying the in- vacuum intrajet character- 
istics, in particular the loss of quantum coherence and the resulting disruption 
of angular ordering for decays inside the medium, the medium-induced radi- 
ation (inelastic energy loss) and elastic scattering of partons on particles in 
the medium (elastic energy loss). Below we shall describe the corresponding 
modifications of the MC procedure for generation of the in-vacuum QCD 
cascades that take these features into account. 

Let us start with the decoherence-induced disruption of angular ordering 
|97j . The violent environment created in ultrarelativistic heavy ion collisions 
acts as a source of random energy-momentum and color with respect to the 
diagrams for the in-vacuum processes. Generally speaking, this kind of ran- 
dom impact ruins the phase tuning lying at the heart of quantum interference 
phenomena. The situations in which random external influences destroy the 
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interference-related effects are very common in solid state physics. A well- 
studied example is provided by the destruction of the weak localization in 
the presence of random external impacts, see e.g. [98| [99] . Thus it is intu- 
itively quite clear that the quantum coherence effects should be broken in the 
violent environment created in heavy ion collisions. This was first explicitly 
mentioned in [M] , in which the origin of decoherence and ensuing disappear- 
ance of angular ordering was related to rescattering of cascading particles 
on scattering centers in the medium. Another argument in support of the 
assumption of angular decoherence is the fact that, at least at the level of 
Vlasov approximation, the characteristic time of color randomization is much 
less than the characteristic scattering time |1UU[ 1101 j . As color randomiza- 
tion alone is sufficient for angular decoherence, one could expect that this 
effect is even stronger than assumed in [91] . The simplest assumption taking 
into account this effect is that angular ordering is broken for all the decays 
that take place inside the medium. The effect is thus directly dependent 
on size of the medium available for the cascade development. In terms of 
notations introduced in Fig. 1, the effect is L-dependent. The corresponding 
modifications of the rapidity distributions are shown in Fig. 2. We see that 
the decoherence effects alone, i.e. before considering medium-induced energy 
loss, lead to substantial softening of the intrajet particle spectrum. 

It is important to notice that, as seen from Fig. 2, the softening of the 
intrajet rapidity distribution sets in at relatively small medium sizes L ~ 
1 fm. 

Let us now turn to the effects of the medium-induced radiation. A rigor- 
ous treatment of this problem is very difficult. Only very recently in |102L 1103] 
an elegant formalism allowing to take into account LPM effect in the MC 
cascades was described. 

The results of modifications of QCD cascade properties were obtained 
using simpler schemes based on various approximations for the effective dif- 
ferential spectrum ( [25] ) or, equivalently, to the modified splitting function 

K tot ( Zj E, Q\ q, L) = K vac (z) + AK(z, Q\ q, L) (46) 

that depends, generally speaking, not only on the energy ratio z but also 
on the energy and virtuality of the parent parton, stopping power of the 
medium, path length in the medium, etc. [104J. 

The simplest version of the medium-modified splitting function mimicking 
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Figure 2: Distribution in rapidity P(y) of final prehadrons: 1) L — Ofm, 
full angular ordering, red, full line; 2) L = 0.5 fm, partial angular ordering, 
green, dashed line; 3) L = 5fm, partial angular ordering, blue, dotted line 



the effects of medium-induced radiation was suggested in [105] : 

K vac (z) -)■ (1 + f me d)K vac (z) 



(47) 



A more advanced definition of the medium-induced splitting function was 
proposed [57] : 

2nQ 2 dl med 



AK(z,E,Q 2 } q } L) 



a s dzdQ 2 



(48) 



where dI med /dzdQ 2 is taken from Eq. (34). 

With the concrete expressions for the modified splitting functions like 
(47), (48) the main ingredient for the description of the cascade, the Sudakov 
formfactor (44), is fully specified. 
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The main quantity addressed in these MC simulations is the jet quench- 
ing ratio Raa(p±)- The agreement with the experimental data is achieved 
by tuning the parameters of the medium like q 2 or f mcd . In most cases the 
agreement can indeed be reached - however, at the expense of very large vari- 
ations in the values of these parameters obtained within different formalisms 
and of the numerous theoretical uncertainties. 

As to the more generic intrajet characteristics, the general expectation 
concerning the effect of the medium-induced interactions are: 

• Softening of the rapidity distributions. 

• Broadening of the distributions in transverse momentum. 

• Increased multiplicity. 

These expectations are indeed confirmed by the studies undertaken with 
the above-described MC generators [9H [95], see also |106| . More involved 
properties of medium-modified cascades such as higher order multiplicity 
moments were also studied in the literature |107| I108j . 

4 The macroscopic approach to the quark- 
gluon medium 

Macroscopic collective properties of a medium may reveal itself in its me- 
chanical motion as a whole (e.g., viscosity) described by hydrodynamics or 
in its response to external color currents (e.g., chromopermittivity) described 
by in-medium QCD. Here we mostly deal with the latter approach discussing 
hydrodynamics at the end very briefly because there exist many extensive 
reviews of it (see, e.g., \W9\ 1110] ). 

Considering the evolution of energy flows due to the elastic scattering, 
bremsstrahlung or synchrotron radiation one follows the fate of an initial 
source (a particle, a parton, a dipole etc) and the medium impact on these 
processes. The evolution of its velocity vector is crucial for conclusions about 
its energy losses. The notion of the coherence length becomes important when 
multiple rescattering is taken into account. 

However, the medium itself can radiate in response to permutations. The 
genuine role of the medium and its collective properties are most clearly 
revealed by its polarization P due to the external current. The macroscopic 
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approach to the description of such collective properties is the most suitable 
one. The linear response of the medium to the electromagnetic field E is 
usually described as 

P = ^E, (49) 

where e is the dielectric permittivity. It is seen that the polarization can be 
quite strong for large values of e. The initial radiation process serves as a 
trigger for the collective response of the medium initiated by the polarization. 
The well known examples are provided by the Cherenkov radiation, the wake 
and the transition radiation. 

What is typical for their description is the (approximate) constancy of 
the particle velocity vector used in the external current. In ultrarelativistic 
processes (7 ^> 1) the relative change of the velocity is much smaller than 
the relative energy loss because they are connected by the formula 

f-(7 2 -D^ (50) 

and the above statement is well supported. The velocity loss can become 
noticeable only for non-relativistic partons. 

In what follows we consider very high energy processes. It is well known 
that the gluons become the main component of the wave functions of the 
colliding hadrons. At LHC, the gg-luminosity is at least by an order of 
magnitude higher than the gg-luminosity. Therefore the in-medium gluo- 
dynamics [11 Ij is considered below. It simplifies the formulae. Quarks can 
be easily included too |112j . 



4.1 Equations of in-medium QCD 

The in-medium equations of gluodynamics differ from the in-vacuum equa- 
tions by taking into account the chromopermittivity of the quark-gluon medium. 
Similar to the dielectric permittivity in electrodynamics it describes the lin- 
ear response of the matter to passing partons. At the leading order in the 
coupling constant the equations are completely analogous to electrodynami- 
cal ones with chromopermittivity just replacing the dielectric permittivity. 
The classical in-vacuum Yang-Mills equations are 

D,F^ = r, (51) 

= d»A v - d u A» - ig[A^, A u ], (52) 
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where = A^T a ; A a (A® = $ a , A a ) are the gauge field (scalar and vector) 
potentials, the color matrices T a satisfy the relation [T a ,T b ] = if a b c T c , = 
dfj, — iglAp, ■), J v {p,S) is a classical source current, and the metric tensor is 
(T=diag(+,-,-,-). 

In the covariant gauge d^A^ = they are written as 

UA v = + ig [ Avf d v A » + (53) 

where □ is the d'Alembertian operator. Thus the leading term in the ex- 
pansion in the coupling constant for the classical gluon field is given by the 
solution of the corresponding abelian problem. 

The chromoelectric and chromomagnetic fields are 

B» = —e>* ij F i: >, (55) 
or, as functions of the gauge potentials in vector notation, 

dA 

E a = -grad$ a - + gf abc A b <S> c , (56) 
B a = curlA a - ^gf abc [A b A c ]. (57) 



The equations of motion (51) in vector form are written as 

divE a - gf abc A b E c = Pa , (58) 
dE 

curlB a - - gf abc (^ b E c + [A b B c ]) = j a . (59) 

The effects due to the medium are accounted for if p a and j a contain 
both external and internal (induced) terms. However, analogously to elec- 
trodynamics, the role of internal currents may be represented in the linear 
response approach by the medium permittivity e if E is replaced by D = eE 
in F^ v , i.e. in Eq. (54). After that is done, only external currents should be 
considered in the right hand sides of the equations. Therefore Eqs (58), (59) 
in vector form are most suitable for generalization to the in-medium case. 
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In terms of potentials the equations of in-medium gluodynamics are cast 
in the form jl 1 II 

d 2 A 1 d 1 

A A a - = ~ 9fabc(-cnr\[A b , A c ] + e— (A 6 $ c ) + -[A fe curlA c ] - 

dA 1 

e*b-^ - e$ b grad$ c - -e// cr)m [A fc [A m A„]] + #e/ cmri $ 6 A m $„),(60) 



<9 2 $ O <9A 9$;, 

A $ a - e°-^ = -7 + ^c(-2A c grad$ 6 + A b ^ - e^$ c ) + 

g 2 famnfnlbA m Ai§ b . (61) 

If one neglects the terms with explicitly shown charge g, one gets the set 
of abelian equations, which differ from electrodynamical equations by the 
color index a only. The most important property of the solutions of these 
equations is that while the in- vacuum (e = 1) equations do not admit any 
radiation processes, for e 7^ 1 there appear solutions of these equations with 
non-zero Poynting vector even in the classical approach. They predict such 
coherent effects as Cherenkov gluons, the wake and transition radiations (see 
section 4.3). This corresponds well to the microscopic description in which 
the matter, at early times after a heavy ion collision, is described in terms of 
the coherent classical field. 

As pointed above, j a is treated now as an external current ascribed to 
partons moving fast relative to the other partons "at rest". It is proportional 



to g as seen from Eq. (67) below. The potentials are also of the order of 



g if one neglects the ^-dependence of e discussed in the section 4.2. Thus 



the terms with explicitly shown g in the right hand sides of Eqs (60), (61) 
would be of the order g 3 . The higher order corrections may be calculated 
(see section 4.5) and quite interesting conclusion about the color rainbow is 
obtained. 



4.2 The chromopermittivity 

The most economical way of description of matter properties is to take into 
account an impact of internal currents in the medium phenomenologically 
by the chromopermittivity e introduced as a constant in Eqs (60), (61). 
Its spatio-temporal dependence may be included similar to electrodynam- 
ics. For the Fourier components of e this would imply the dependence on 
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the frequency u and the wave vector k. It is well known in electrodynam- 
ics that the magnetic permeability is then automatically taken into account. 
Therefore it is necessary to consider the permittivity tensor and its color 
dependence. 

These complications can become important but even in electrodynam- 
ics they are often ignored. The medium is often assumed to be isotropic 
and homogeneous. Methods of linear response theory are applied. The in- 
dependence is commonly accounted for because the dielectric permittivity 
of ordinary substances depends on u. Moreover, it has an imaginary part 
determining the absorption. The energy behavior of e determines the scales 
specific for a particular substance. 

All these features are usually found experimentally using the relation 
between the dielectric permittivity and the refraction index n 

e = n 2 (62) 

For example, both Ree and Ime for water have been measured in a wide 
energy interval ranging over 20 orders of magnitude |113] . The real part Re e 
is approximately constant in the visible light region (-^/e ~ 1.34), increases 
at low uj (up to e ~ 80) and becomes smaller than 1 at high energies, tending 
to 1 at frequencies exceeding the Langmuir (plasma) frequency as 



to 2 

Ree = 1 - -4. (63) 
ur 



The imaginary part (Ime) responsible for the absorption is very small for 
visible light but drastically increases in nearby regions both at low and high 
frequencies (in this way the Nature saves our eyes!). 

Thus there exist at least two characteristic scales for Ree: just below the 
optical region and the high energy scale determined by the Langmuir fre- 
quency related to the proper plasma oscillations. The borders of the optical 
region define the scales for Ime. Theoretically this behavior is ascribed to 
various collective excitations in water relevant to its response to radiation 
with different frequencies. Among them the resonance excitations are quite 
prominent (see, e.g., |114j ). Even in electrodynamics, the quantitative the- 



ory of this behavior is still lacking, however. Moreover the formula (63) is 
purely electrodynamical one and does not take into account hadronic pro- 
cesses at extremely high energies where the hadronic components of photon 
wave functions are important. 
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Then, what can we say about the chromopermittivity? 

The chromopermittivity can be calculated from the polarization operator. 
Up to now, the attempts to calculate the chromopermittivity from first prin- 
ciples are not very convincing. The corresponding dispersion branches have 
been computed in the lowest order perturbation theory 1116] \117\ 1118] . 
The properties of collective excitations have been studied in the framework of 
the thermal field theories [TEHJ 11201 H2D Q221 CCS] . In these works the simplest 
perturbative results representing e are generalized by summing some types of 
diagrams. Only the time-like dispersion relations were obtained in this way. 
The space-like dispersion relation was obtained in |124] in the framework of 
the method using an ad hoc assumption about a special role of scalar reso- 
nances. Some studies were performed in |112j . All of these works refer to the 
weakly coupled medium. The phenomenological model built upon some in- 
tuitive guesses (see, e.g., |125j ) were also studied. The physical consequences 
of introducing the chromopermittivity could be mimiced in microscopic ap- 
proach considering different non-zero masses and coupling constants of in- 
vacuumand in-medium gluons (see, e.g., [72| l?3] 174] fl2T) ] IT271 fT28l fT29] 1X30] ) . 
That allows to vary the dispersion relations. However, no convincing prin- 
ciples of their choice exist and the number of adjustable energy dependent 
parameters increases even though they are related with more fundamental 
microscopic notions. 

Experimental data surely favor strong response of the medium to de- 
posited energy. It implies that the chromopermittivity must be defined in 
terms of non-perturbative in-medium gluon field correlation functions that 
are as yet unknown. Therefore, it makes sense to consider the chromoper- 



mittivity as a complicated function of g and, correspondingly, the Eqs (60) 



(61 ) as strongly non-linear in g even at the classical level because e ^ 1 takes 
into account some quantum and non-abelian effects. 

In view of this situation, we prefer to use the general formulae of the 
scattering theory |113[ 1131 j to estimate the chromopermittivity. It can be 
expressed |113] 11311 1132] through the real part of the forward scattering 
amplitude ReF (w) for the corresponding quantum^] 

ReAe = ReeM - 1 = *W-(") = (M) 

CO 2 u 

with 

ImFo(w) = ^-a(u). (65) 
4tt 

11 In electrodynamics these quanta are photons, in QCD they are gluons. 
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Here u denotes the energy, N s the density of scattering centers, u(co) the 
cross section and p(oj) the ratio of real to imaginary parts of the forward 
scattering amplitude F (u). The sign of Ae coincides with the sign of p{oS). 

Unfortunately, we are unable to calculate directly in QCD these charac- 
teristics of gluons and have to rely on analogies and our knowledge of the 
properties of hadron interactions. The experimental facts we get for this 
medium are brought about only by particles registered at the final stage. We 
hope that some of their features are also relevant for gluons as the carriers of 
strong interactions. Those are the resonant behavior of amplitudes at rather 
low energies and the positive real part of the forward scattering amplitudes at 
very high energies for hadron-hadron (pp, Kp, up) and photon-hadron ('jp) 
processes as measured from the interference of the Coulomb and hadronic 
parts of the amplitudes. The latter feature is directly connected with the 
specific property of hadronic processes, the energy increase of the total cross 
sections, by the dispersion relations for forward scattering amplitudes. 

Any Breit-Wigner amplitude has a positive (negative) real part below 
(above) its peak. The combined effect of all collective resonance excitations in 
the medium may lead to positive Ae in a wide energy interval at low energies. 
This is the necessary condition for Cherenkov effects. As explained in the 
paragraph 4.3.1, the double-humped events at RHIC and the asymmetry of 
mass shapes of resonances passing through the quark-gluon medium observed 
at SPS may be explained in this way. 

At very high energies above the threshold ojth of the positiveness of real 



parts of the amplitudes in (64) one could try to use a simple model with 



Ree=(l + ^)e(w-a; tft ). (66) 



The model (66) satisfies the Kramers-Kronig requirement for Ree to be an 



even function and approaches 1 at infinity. If experimental data on Re_F (u;) 



for hadrons are inserted in Eq. (64) one gets the shape of Ree very close to 
Eq. (66) (see Fig. 1 in [1331 1134] ). The cosmic ray events with the ringlike 
(conical) structure may be then explained. 

At intermediate energies in between these two intervals the similar argu- 
ments would favor negative values of Ae. Thus, no Cherenkov gluons would 
be emitted with such energies. The transition gluon radiation (section 4.3.3) 
is not forbidden, however. The wake effects (section 4.3.2) can be shown to 
be small (see the paragraph 4.3.2). 
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It is tempting to conclude that energy scales determining the behavior of 
the real part of the chromopermittivity are not directly related either to the 
temperature T or to the QCD scale Aqcd of the order of hundreds MeV but 
extend to higher values defined by the end of the resonance region (about 
several GeV) and by the threshold value u t h connected with energy increase 
of hadronic total cross sections (about tens of GeV). The microscopic nature 
of their origin is yet unclear. 

There is another important problem of the definition of the chromoper- 



mittivity. The in-medium equations (60), (61) are not Lorentz-covariant. In 



macroscopic electrodynamics the rest system of the matter is well defined 
and the dielectric permittivity is considered just there. For collisions of two 
nuclei (or hadrons) the definition of e depends on the experiment geometry. 
Therefore one needs to consider the Lorentz-covariant tensor of the chromop- 
ermittivity. The spectra of classical Cherenkov radiation in moving media 
were calculated in |135j and with quantum corrections in |136j . 



The notions of dilute and dense systems (i.e. the value of N s in (64)) de- 
pend on values of x and pt at which the wave functions of colliding nuclei are 
probed. Due to Lorentz transformations, partons moving in different direc- 
tions with different energies can "feel" different states of matter in the same 
collision of two nuclei because of the dispersive dependence of the chromoper- 
mittivity. We discuss this when dealing with different experimental data and 
with the problem of instabilities of the quark-gluon medium (section 4.4). 

4.3 Classical polarization effects in the quark-gluon 
medium and its chromo dynamical properties 

In view of the similarity between classical QCD and QED equations it would 
be of no surprise to observe some effects in nucleus collisions reminding those 
in electrodynamics. Actually, this idea was first promoted in |133l 1134] re- 
lying on the similarity of quarks to electrons and gluons to photons. The 
emission of Cherenkov gluons analogous to Cherenkov photons was predicted 
and used for interpretation of a cosmic ray event just observed [137J. These 
events are discussed in the paragraph 4.3. l.b. The strong support this idea 
got, however, much later from RHIC data on nucleus-nucleus collisions (see 
the paragraph 4. 3.1. a). 

Before delving in experimental data let us describe the classical solutions 



of Eqs (60), (61). As usually, the current with constant velocity v along the 
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z-axis is considered: 



j(r,t) = vp(r,t) = 47rgv<5(r — vt). (67) 

The color index is omitted because it points in a fixed direction in color group 
space. The classical lowest order solution of in-medium gluodynamics can be 
cast in the form I138j 

$ (i) (r t) = 2g e(vt-z-r ± Vev 2 -l) 
e y/(vt-z) 2 -r 2 ± {ev 2 -iy 

and 

A«(r,£) =ev$( 1 )(r,t), (69) 

where the superscript (1) indicates the solutions of order 0(g) (not to count 
for g-dependence of e), rj_ = ^ x 2 + y 2 is the cylindrical coordinate; 2 is the 
symmetry axis. These formulae describe both the space-time profile of the 
energy-momentum deposition and the dynamical response of the medium 
(the value of e). 

This solution describes the cone-like emission of Cherenkov gluons at the 
typical angle 

cos# = ^. (70) 

It is constant for constant e > 1. 

The expression for the intensity of the radiation is given by the Tamm- 
Frank formula [159] (up to the Casimir factors Cr) 



dE 



47rasC R [ udu(l - -^—)&(v 2 e(io) - 1). (71) 
J v z e(u) 

For absorbing media e acquires the imaginary part. The sharp front edge 



of the shock wave (68) is smoothed. The angular distribution of Cherenkov 



radiation widens. The S- function at the angle (70) is replaced by the a'la 
Breit-Wigner angular shape |140[ 11411 1142] with the maximum at the same 
angle (for small imaginary part of e) and the width proportional to the imag- 



inary part (see Eq. (78) below). The special role of 1/e term in (68) which 
describes the wake left behind the parton is discussed in section 4.3.2. 

Let us stress that this classical effect has collective non-perturbative origin 



even though as enters seemingly linearly in Eq. (71). The chromopermit 



tivity e takes into account the non-perturbative terms responsible for the 
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collective medium response. Only the emission of a primary gluon which 
triggers this response is treated perturbatively, hence the factor of as in Eq. 



(71 ), but its effect is strongly enhanced and modified by the medium as seen 



from the term in brackets in the integral. 



4.3.1 Cherenkov gluons 

a. The double-humped structure at RHIC 

The scalar and vector potentials in the momentum space can be cast in 
the form 

5(u - kvQv 2 ? (?2) 



= 2ngQ a 



uj 2 e(ev 2 ( 2 
A (i) = $(i) 
( = cos 9, 



(73) 
(74) 



u, k are the energy and momentum, 9 is the polar angle. Let us note again 
that both the Cherenkov term (ev 2 ( 2 — 1) _1 and the wake factor e -1 are 
present. 

The energy loss dW per the length dz is determined by the formula 



dW 

dz 



-gE z . 



In the lowest order 



d 4 k 



[uA^\^uj) -A; 2 $ (1) (k,u;)]e i(kv - 



-Lj)t 



(75) 



(76) 



Inserting (72), (73) in (76), ( |75j ) one gets (see also |143j ) 
dW^ g 2 uj f v 2 h-( 2 ) l 



dzdQdu) 2ir 2 v 2 ( 



1 - e t v 2 ( 



(77) 



The first term in the brackets corresponds to the transverse gluon Cherenkov 
radiation (index t at e t ) and the second term to the radiation due to the lon- 
gitudinal wake (index I at e{). The transverse and longitudinal components 
of the chromopermittivity tensor are explicitly indicated here even though 
they are equal in any homogeneous medium. 
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Similar to electrodynamics |14Uj . one easily gets from (77) the energy- 
angular spectrum of emitted gluons |141[ l!42j per the unit length 



where 



x 



c 2 , 



dVtduj 



%o — tit/ 



a s 



2tt 5 



I |2 2 



x)Y t 



+ 



(x-x ) 2 + (r t ) 2 /4 x 



r, = 2e 



2.;/ 



|2„.2 



(78) 



(79) 



The real (ei) and imaginary (e 2 ) parts of e are taken into account. The 
angle 6 is the polar angle if the away-side jet axis would be chosen as z- 
axis. It is clearly seen from Eq. (78) that the transverse and longitudinal 



parts of the chromopermittivity are responsible for the distinctly different 
effects. The ringlike Cherenkov structure (conical emission) around this axis 



is clearly exhibited in the first term of ( 78 ) . The second term defined by the 



longitudinal part of e is in charge of the wake radiation described in section 
4.3.2. 

Let us consider the first term. The angle 9 is related to the laboratory 
polar (Ol) and azimuthal (</>/,) angles in RHIC experiments through 



x = cos 6 



sin 2 Q L cos 2 ■ 



(80) 



Integrating ( 78 ) over $l one gets |141j the final formula to compare with the 
two- hump structure of azimuthal (4>l) correlations observed at RHIC. The 
formula is quite lengthy and therefore not reproduced here. It is seen already 



from (78) and (80) that this projection of the two-dimensional ring on its 
diameter is symmetrical about 4>l = tc and exhibits humps whose positions 
are mostly determined by e\ and widths determined by e 2 . The first term in 



(78) clearly displays the Breit-Wigner angular hump which replaces the 5- 



functional angular dependence characteristic for real e (easily obtained from 



(J78J) at r t -»■ 0). 

Actually, it is a matter of experiment to demonstrate whether Cherenkov 
gluons can be observed in some energy regions. No definite knowledge about 
the chromopermittivity and its excess over 1 except the above conclusions is 
available. Moreover, the absorption may be strong enough and hadroniza- 
tion (confinement) may change the final characteristics. The very first ex- 
perimental indication in 1979 |137j which favored this effect dealt with a 
single cosmic ray event (also some others less pronounced were observed in 
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cosmic rays and at lower energies in emulsion experiments at accelerators 

[HI Eg EH [T13 

That is why it was a surprise to find in 2004 and confirm later [158, 
11551 11601 fTBTl H621 IT631 fTM flfel IT661 IT6T1 [168] the so-called double-humped 
events in central nucleus-nucleus collisions at RHIC. In fact, the obtained 
distributions closely remind those published by Cherenkov in 1937 when he 
used the projection of famous rings on their diameter and got the double- 
humped structure (e.g., see [169J). Let us stress that this effect reflects 
an intrinsic property of the medium phenomenologically described by its 
permittivity e. 

Experimentalists at RHIC used a quite elaborate method to learn about 
a parton passing through the quark-gluon medium. Imagine the large angle 
scattering of two partons belonging to the colliding nuclei which happens 
near the surface of the medium. The triggers placed at the angle about n/2 
to the collision axis would register a normal jet created by a parton pass- 
ing through the thin layer. Another scattered parton moves through the 
quark-gluon medium. Both partons produce jets. The two-particle correla- 
tions were studied. Those particles which belonged to the same jet formed a 
peak at low A0. It was expected to get a peak in the opposite direction at 
A0 m 7r due to the correlation between particles from the trigger (near-side) 
and away-side jets. This was really found in pp-collisions but not in cen- 
tral nuclei collisions. It implied that the away-side jet is strongly modified 
when passing through the quark-gluon medium. For pions produced in some 
energy intervals it could not be simply attributed to jet quenching because 
the two-hump structure appeared. Namely this phenomenon was related to 
Cherenkov gluons. 

In |141] a fit of this experimental data according to the formulae de- 
scribed above was performed and the values of real and imaginary parts of 
the chromopermittivity were obtained. The Monte Carlo program taking into 
account the specifics of collisions of initial partons and hadronization at the 
final stage (described by the transverse momentum spread parameter Aj_) 
was developed. The fit to experimental data (with elliptic flow subtracted) 
is shown in Fig. 3. 

In this way it is possible to determine the values of all the three adjustable 
parameters (see Table I). The constant values of were used because the 
interval of the transferred momenta (frequencies) studied in experiment is 
quite narrow. Moreover, this assumption is directly supported by experiment 
(for recent results see |162[ \1Q3\ 11641 1170j ) according to which the positions 
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Figure 3: Associated azimuthal correlations at STAR: circles - experiment, 
triangles- theory. 



of peaks do not depend on transverse momenta of both the trigger (up to 
pip 9 = 6 GeV) and associated particles (up to p9f soc = 4 GeV). That implies 
constant e\ according to (78), (79). 



Table 1 



Experiment 


^max 


ei 


£2 


A_l, GeV 


New data 


STAR 


1.04 rad 


5.4 


0.7 


0.7 


6> max rad 


PHENIX 


1.27 rad 


9.0 


2.0 


1.1 


ei « 6; e 2 w 0.8 



The initial results differed for STAR and PHENIX because of disagree- 
ment in peaks positions (mostly determined by ej), tt ± 1.04 and ir ± 1.27, 
correspondingly, and the widths (mostly determined by e 2 ). The new data 
|170[ 1171] give 7r ± 1.1 so that e\ ~ 6 and e 2 ~ 0.8 would be good estimates 
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as can be guessed from the range of values in Table 1. The main conclusion 
about the large values of e\ indicating on the non-gaseous matter with large 



N s in Eq. (64) is supported in any case. The imaginary part of the chromop- 
ermittivity is comparatively small relative to the real part. Low attenuation 
helps observe this effect. 



The experimental data show that the often citec 12 conclusion of [124 J 
about smaller angles for larger px is wrong. The model of a single scalar 
resonance used there for the derivation of the dispersion equation defining 
the energy dependence of the chromopermittivity is oversimplified. This 
defect can be cured by the set of overlapping resonances from PDG which 
leads to approximately constant e in the low-energy region cut off by the 



Hagedorn temperature if the formula (64) is used. This set of resonances 
can be used to model the collective excitations of the quark-gluon medium 
analogously to excited levels of atoms in electrodynamics |114j . The energy 



loss per 1 Fm can be estimated according to Eq. ( 71 ) as about Cr GeV 
for the end of the resonance region near 4 GeV. It depends quadratically 
on this parameter that induces some flexibility in the estimate which shows, 
nevertheless, that they are rather large. 

For three-particle correlations |171|. \172\ 1173] the peak, according to the 
initial data, was at large angles. However, the latest data |171j agree with the 
above-described estimates. These correlations show also that the hypothesis 
about the deflected jets |174j is not valid because in experiment there are 
off-diagonal peaks which can not appear according to this hypothesis. The 
presence of such peaks is an evidence for the conical emission pattern. 

In- medium modifications of QCD bremsstrahlung calculations (see, e.g., 
|175j ) predicted some broadening of the away-side peak which was, however, 
not enough to get the double-humped structure after subtracting the collec- 
tive flow V2 effects. 

In principle, the two-hump structure may arise not only due to Cherenkov 
gluons but also as the sonic Mach cone |176[ \177\ 1178] . In distinction to 
gluodynamics, it is ascribed to the longitudinal excitations in plasma |178] 
125J. The model with an hypothesized space-like dispersion relation 



uj = ^/u 2 k 2 + uol (81) 

was considered in |125j . Here ui p denotes the plasma frequency, u is the speed 
of a plasmon. It describes a Mach shock wave for a supersonical source. The 
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E.g., "the Cherenkov angle decreases with the momentum of the radiated gluon" |38j . 
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similar structures appeared in the models using the AdS/CFT conjecture 

pzaiiism. 

In the 3D hydrodynamics with big local density fluctuations which at the 



initial moment have the form of the tubes _ there arise both the ridge and 
the double-humped structures |181[ 11821 1183] . They appear only if the com- 
bination of these special non-smooth initial conditions with hydrodynamical 
evolution of the hot and dense medium is used so that the tubes induce an 
appearance of the ridges in the hydrodynamical flow. At the same time, 
as was shown in |184j . in hydrodynamics with somewhat different and less 
special initial conditions the particle yield from the transverse wake moving 
opposite to the trigger jet largely overwhelms the weak Mach cone signal for 
relativistic particles ^ The jet of particles in this direction should appear 



(the strong away-side jet) that is not observed in experiment. As shown in 
section 4.3.2, this effect does not take place in chromodynamics where the 



longitudinal wake (the second term in Eq. (78)) results in emission perpen- 
dicular to the parton trail. The hydrodynamical Mach cone interpretation 
was not yet used for quantitative fits of experimental data aimed at getting 
information on the medium properties. 

The charge's orientation in color space is fixed in classical equations of 
gluodynamics. There is no color dependence in their solutions beside trivial 
Casimir factors which just determine the overall normalization. Thus the 
statement that "their absence ... in "conical flow" suggests sound radiation 
rather than the gluonic one" |185j is absolutely unjustified. The functional 
color dependence may appear in gluodynamics as an effect at higher orders 
(see section 4.5). 

There were several attempts to ascribe the double-humped structure 
around away-side jets to the large angle emission of shower gluons |174[ I185j 
or to their deflection due to random multiple scattering |186j or coupling to 
collective flows |187l 1188] but they ask for additional ad hoc parameters with 
values that are not compatible with experiment |162[ I163[ 1164] . The collec- 
tive medium behavior as the origin of humps is clearly seen from properties of 
hadrons inside the shoulder regions (pt- and centrality-independent shapes, 
mean p?, composition) being similar to those observed for inclusive processes 



13 Does it not remind the field tubes in Glasma described above? 

14 This conclusion differs from those in [181, 182, 1183] . where no away-side signal at 
A<f> — 7r is seen. Probably, this is due to the choice of a single high density spot in 
initial conditions of the hydrodynamical model while two partons are necessarily created 
in parton treatment. 



50 



|162[ 11631 1164] . Their energy dependence suggests stronger medium effects 
at higher energies, as expected in [1621 11631 1164j . 

The method of wavelet analysis was successfully applied to distinguish 
between jets and conical structures in three-dimensional plots on the event- 
by-event basis [1551 IT5D] . 

It was assumed everywhere above that the direction of gluons defines the 
direction of final pions. The new mechanism of direct emission of Cherenkov 
mesons by heavy quarks was considered in the effective theory with account 
of the gauge/gravity duality [127J. The Cherenkov angle is the same, but 
it is determined now by the relation of masses i.e. by the kinematics alone 
without any reference to medium properties (the chromopermittivity). 

b. The asymmetry of in-medium resonances (SPS etc) 

The necessary condition for Cherenkov effects to be observable within 
some energy interval is e(uS) > 1 for u belonging to this interval. Accord- 



ing to Eq. (64) this implies that ReF (w) > 0. This requirement is ful- 
filled within the low-energy (left) wings of resonances described by the Breit- 
Wigner formula. Therefore one could expect that the collective excitations of 
the quark-gluon medium leading to Cherenkov gluons contribute in these en- 
ergy intervals in addition to the traditional resonance effects. Herefrom one 
gets the general prediction that the shape of any resonance formed in nucleus 
collisions must become asymmetrical with some excess within its left wing 
compared to usual Breit-Wigner shape. Since the probability of Cherenkov 
radiation is proportional to Ae = e — 1 this asymmetry must be proportional 
to p(oj). Experimentally one of the best ways to observe it is to measure, 
say, the dilepton mass distribution of the corresponding modes of resonance 
decays. Then this distribution must have the shape |191] 

dNu A ( m 2 -M 2 ^, A 

'l + w r r — 9(m r -M) (82) 



dM (m 2 - M 2 ) 2 + M 2 r 2 V M 2 

Here M is the dilepton mass and m T is the resonance mass. The first term 
corresponds to the ordinary Breit-Wigner cross section with a constant nor- 
malization factor A. According to the optical theorem it is proportional to 
the imaginary part of the forward scattering amplitude. The second term de- 
scribes the coherent Cherenkov response of the medium proportional to the 
real part of the amplitude. Its contribution relative to ordinary processes is 
described by the only adjustable parameter w r for a given resonance r which 
must be found from comparison with experimental data . 
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This formula has been compared with the most accurate experimental 
results |192j on the dilepton spectra of the p-resonance produced in nuclei 
collisions at SPS (see Fig. 4). The asymmetry of dilepton mass spectrum 
is clearly seen. The corresponding parameter w p was found to be equal to 
0.19. This confirms the discussed above conclusion that the contribution of 
Cherenkov effects is not negligible and stresses the fact that resonance regions 
may be responsible for it at rather low energies. This feature is well known 
in electrodynamics (see, e.g., Fig. 31-5 in |114j ) where the atoms behaving 
as oscillators emit as Breit-Wigner resonances when get excited. 




Figure 4: Modified spectrum of dileptons in semi-central collisions In(158 A 
GeV)-In for NA60 (points) compared to the p- meson peak in the medium 
with additional Cherenkov contribution (dashed line). 

Nowadays there exist numerous data on other resonances produced in 
nuclei collisions with qualitative indications on asymmetry of their shapes 
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with excess in low- mass wings |193[ I194[ I195[ 11961 1197j . This agrees with the 
universal prediction of Eq. (82). 



c. Cosmic ray events 

The very first events with the ring-like structure were observed in non- 
trigger cosmic ray experiments. Namely they initiated the idea about Cherenkov 
gluons [133J. Two peaks more densely populated by particles than their sur- 
roundings were noticed in the cosmic ray event |137] initiated by a primary 
with energy about 10 16 eV close to LHC energies. 



n n r~n 




-/ f Z J 4- 5 Lnr t - 

Figure 5: Distribution in the number of produced particles at different dis- 
tances from the collision axis r in the stratospheric event at 10 16 eV [137] 
has two pronounced peaks. Correspondingly, the distribution over pseudora- 
pidity also has two such peaks. 



The distribution in the number of produced particles is shown in Fig. 
|5j where it is plotted as a function of the distance from the collision axis 
proportional to the polar angle 9. It clearly shows two maxima. This event 
has been registered in the detector with nuclear and X-ray emulsions during 
the balloon flight at the altitude about 30 km. At the very beginning the idea 
about clusters (fireballs, clans, jets etc) was proposed but it failed to explain 
the data. Approximately at the same time the similar event with two peaks 
was observed at 10 13 eV |144j . Some events with one peak (due to the limited 
acceptance of the installation?) were shown even earlier |145l 11461 1147] . 
When the two-dimensional distribution of particles was considered in the 
azimuthal plane (called as the target diagram in cosmic ray experiments), 
this event revealed two (forward and backward in c.m.s.) densely populated 
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ring-like regions within two narrow intervals of polar angles (corresponding 
to peaks in Fig. [5]) but widely distributed in azimuthal angles for each 
of them. Therefore such events were termed ring-like events. The peaks 
were interpreted as effects due to Cherenkov gluons emitted by the forward 
and backward moving initial high energy partons. In cosmic ray events the 
energies of partons are high and the effect can be related to the high-energy 
behavior of KeFo(E). 

It is important to notice that the energies of particles produced within 
the humps (rings) in RHIC data are rather low on the scale of initial energies. 
Even the energies of jets-parents in the double-humped events are lower than 
5 GeV. Therefore the whole effect at RHIC can be related to the low-energy 
region. 

One of the most intriguing problems is to understand properly the fact 
that the RHIC and cosmic ray data were fitted with very different values 
of the chromopermittivity, close to 6 and 1 correspondingly. This could be 
interpreted as due to the difference in values of x (the parton share of energy) 
and Q 2 (the transverse momenta squared) and the Lorentz-transformation 
to different rest systems. It is well known that the region of large x and Q 2 
corresponds to the dilute partonic system. At low x and Q 2 the density of 
partons is much higher. Thus the effective values of N s in (64) are different. 

As follows from the estimate based on Eq. ( 64 ) |198l 1199] , the density 
of scattering centers at RHIC conditions is very high, about ten partons per 
proton volume. In this case one deals with rather low x and Q 2 . Therefore, 
one should expect the large density of partons in this region and large Ae. It 
is interesting to note that the two-hump structure disappears in RHIC data 
at higher px where the parton density gets much smaller. It corresponds to 
smaller e and 6, i.e. humps merge in the main away-side peak. 

In the cosmic ray event the observed effect is related to the leading partons 
with large x. Also, the experimentalists pointed out that the transverse 
momenta in this event are quite large. In this region one would expect for 
low parton density and very small Ae. In distinction to a single parton 
traversing the medium at RHIC energies, the collision of forward moving 
partons at the LHC may be considered as the collision of the two relatively 
dilute bunches of partons. Nevertheless, the macroscopic approach is valid 
there as well. 

Thus the same medium can probably be considered as a liquid or a gas 
depending on the parton energy and transferred momenta. This statement 
can be experimentally verified by using triggers positioned at different angles 
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to the collision axk 15 and considering different transverse momenta. In that 
way, the hadronic Cherenkov effect can be used as a tool to scan (l/x,Q 2 )- 
plane and plot on it the parton densities corresponding to its different regions. 

4.3.2 The wake 

In central collisions the interaction region is symmetrical and the change of 
the trigger position might lead to some effect only if the chromopermittivity 
strongly depends on w. It is not likely to be the case for narrow energy 
intervals. However, in mid-central nuclear collisions the interaction region 
itself becomes asymmetrical reminding the rugby ball. Therefore new effects 
can appear when the angular position of the trigger is changed. It has been 
found (see Figs 2 and 3 in |171j ) that the similar two- hump structure appears 
but with some additional contribution noticeable at a particular orientation 
of the trigger particle just in between the in-plane and out-of-plane direc- 
It was stated in |171j that "at present, it is unclear whether this 



tions 
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merely reflects a geometry dependent shift in the away-side peaks or perhaps 
an additional contribution at A<p = 7r/2". 

We have provided arguments [142J that this is an additional contribution 
in the vicinity of A0 = 7r/2 known in electrodynamics as the wake effect 
(see, e.g., pOOj). At the same time, it is strongly influenced by the special 
geometry of mid-central collisions in such a way that the wake radiation at 
about 7r/4-orientation of the trigger particle can escape the overlap region 
much easier than wake gluons for in- and out-of-plane orientations of the 
trigger particle. 



According to Eq. (78) the second term corresponds to the radiation 
induced by the longitudinal oscillations in the wake left behind the away- 
side parton. It is emitted mostly in the direction perpendicular to the wake. 
That reminds the dipole radiation. Namely the color charge-changing regions 
of the dipole-like structure behind the parton are seen in the Monte Carlo 
simulations of the wake [20111202] trailing the color parton. The enhanced and 
depleted charge density regions alternate. The 1/ A/x-singularity is however 
stronger than the sin 2 6 dipole distributions. It must be somehow saturated 
in the vicinity of x = by the spatial dispersion of e. Therefore one may use 
this term for qualitative conclusions in the region not very close to x = 0. 



15 An example is described in the next section 4.3.2. 

16 This plane is defined by the collision axis and the impact parameter. 
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The estimates [142] of the relative importance of the two terms in Eq. 
(78) show that the wake radiation is negligible at the peak of the Cherenkov 



humps at Xq. However, they become comparable at 

*• " 2^Tl < 83 > 

that corresponds to the angles tt — A<pL « 1.43 rad. The wake radiation 
overwhelms the Cherenkov contribution at A(j)£ < n — 1.43 rad where the 
latter decreases. Therefore the shift of their combined maximum to 7r — 
A(/)l ~ 1.3 rad noticed in |171] and shown in Fig. 2 in this paper seems quite 
plausible if explained as the effect of the wake. 

Let us note the difference between Cherenkov and wake radiations. Some- 
times the name of wake is attributed to both of them |125[ 1178] . However, 
it can be shown |142] |200j that the trail behind the parton (the wake) in- 
duced by the chromopermittivity corresponds just to the second term in ( J78~| ). 
Moreover, the wake radiation exists only for non-zero imaginary part of 6/ 
while the Cherenkov radiation also exists and gets the traditional ^-functional 
angular shape for real e. The effect of transverse wake in hydrodynamics is 
completely different. It produces a strong peak in the away-side direction at 
A<pL = Ti" [184] which is not observed in experiment. 



4.3.3 Transition radiation 

In 1945 V.L. Ginzburg and I.M. Frank proposed the idea about the transition 
radiation |203] . Any electric charge passing from some medium to another 
one or moving with (almost) constant velocity in inhomogeneous medium 
with variable dielectric permittivity induces its polarization and radiation of 
photons. This property has been used, e.g., to build the transition radiation 
tracker in ATLAS (A Toroidal LHC Apparatus) detector at the Large Hadron 
Collider (LHC). The transition radiation is determined by restructuring of 
the electromagnetic field surrounding the electric charge. 

During hadron (nucleus) collisions the partons intersect the surface of 
the partner particle (nucleus). This transition should induce abrupt changes 
of their color fields and, consequently, radiation of gluons. The classical 
in-medium QCD equations coincide in their general structure with electro- 
dynamical equations. Since this approach happens to be successful with 
Cherenkov gluons and wakes it is tempting to apply the Ginzburg- Frank for- 
mula to estimate the number of gluons emitted per unit length within some 
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frequency and angular intervals 
dN a s Je^sm 2 26 



dudVL ixuj 

(e 6 - e a )(l - e a + y/ e h - e a sin 2 6) 



[1 - e a cos 2 6>)(l + \fe b - e a sin 2 0) (e 6 cos 8 + A/e a e 6 - e^sin 2 ' 



(84) 



where e a ^ are the chromopermittivities of the media a, b as felt by relativistic 
(v ~ c) partons traversing the surface when passing from b to a. In equation 



(84) we have omitted the color Casimir factors. 

There is important difference between Cherenkov and transition radia- 
tions. Cherenkov radiation can be observed only if the particle moves in the 



medium with e > 1. Its intensity is proportional to e — 1 according to (71 ). 
Transition radiation is proportional to (e^ — e a ) 2 (see (84)) and appears at 



any change of e. In electrodynamics it is studied mostly in the frequency 



interval where the formula (63) is applicable i.e. where Ae < 0. Namely this 
behavior of the dielectric permittivity is responsible for specific features of 
the transition radiation. In the case of the nuclear medium it is unknown 
whether the analogous formula may be used in any energy region. For nar- 
row energy intervals it is reasonable to assume Ae to be constant. At very 



high energies the dependence (66) may be used. Both these assumptions find 
some support in RHIC and cosmic ray data as discussed above. That is why 
we consider [204J first these two possibilities. 

The away-side parton at RHIC is created inside the quark-gluon medium, 
moves through it and escapes in the vacuum where it hadronizes. Accord- 



ing to the above estimates e& ~ 6 ^> e a — 1. Using (84) one gets the 
bremsstrahlung-like spectrum with both infrared and collinear divergencies 
in asymptotics 

dusdO 

dN oc — — . (85) 

L) 

At finite energies the angular divergence is replaced by the maximum at 
9 Ki y 1 . Anyway, this radiation is collimated near the initial direction of 
the away-side parton and is hard to detect because of the strong background. 
The similar conclusions are obtained |205] for the radiation by a gluon moving 
in a stochastic medium. It is interesting to note that for e& 3> e a the intensity 



of the transition radiation (85) does not depend on e b - 



According to (66) for the forward moving partons with extremely high 
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energies the deviation of chromopermittivity from unity is very small. Then 

dN X ^Q( U - Uth ). (86) 

The collinear divergence remains the same but the energy spectrum gets a 
strong peak near u t h- It implies that the resonance-like almost monochro- 
matic subjets with masses close to u t h will be produced at LHC energies. 
These "resonances" are mostly created by gluons and, therefore, are neutral. 
The most favored channel would be to observe them as peaks in mass spectra 
of forward moving -pairs near u t h- Unfortunately, there exists no clear 
prescription how to translate the threshold energy in hadronic reactions to 
gluons. If observed, peaks locations and their disappearence at u — u t h due 
to Q(u — Uth) would provide us information about that. 

One may hope that the confinement and hadronization do not spoil these 
conclusions. This is supported by observation of effects due to Cherenkov 
gluons. However, the transition radiation might be stronger influenced by 
confinement near the QGP surface than collective excitations inside the vol- 
ume. From one side, it shortens the radiation length but, from another side, 
the abrupt variation of color fields surrounding partons when they come out 
of the quark-gluon plasma and hadronize may enlarge the transition radiation 
because of strong confining forces. At the same time the transition radia- 
tion may appear in much wider energy intervals than Cherenkov radiation 
(in particular, at u < Uth)- However, neither experimental nor theoretical 
arguments about the chromopermittivity in these regions exist nowadays. It 
would not be surprising if we find its effects outside the regions where the 
double-humped events were observed. That may somewhat smoothen down 
the structures near uj th . 

Another way to account for difference of vacuum and medium chromoper- 
mittivities is to consider different masses of gluons and values of the running 
coupling in these two enviroments. It was done in |126[ I128[ I129j where the 
transition radiation of heavy quarks was treated in microscopic approach. 
It was shown that it can be as important as the induced gluon radiation. 
However, the conclusions are still indefinite because of strong cancellation 
effects and assumptions about the behavior of the coupling constant. 

One should also point out important problem of the finite size L of 
hadronic objects. It always induces (see, e.g., [126] I130[ 11341 1206] ) the factors 
like 1 — cosL/Lq in emission probabilities. Estimates of the scale length L 
are still not very definite. This problem is closely connected with decisions 
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about the radiation (formation) lengths and/or about the points where the 
partons were created. 



4.4 Instabilities at high energies 



The a'la tachyonic model (66) proposed for the chromopermittivity at high 



energies is inevitably related with instabilities of the quark-gluon medium. 



The dielectric permittivity of the macroscopic matter (e.g., as given by (63)) 
is usually considered in its rest system which is well defined. For collisions 
of two nuclei (or hadrons) such a system requires special definition (see |1114 
207J). In particular, for fast forward moving partons the spectators (the 
medium) are formed by the partons of another (target) nucleus at rest. Thus 
we consider a problem of the system of the quark-gluon medium impinged 
by a bunch of fast partons similar to the problem of plasma physics, namely, 
that of the interaction of the electron bunch with plasma [208, 209l l21Uj . This 
differs from RHIC conditions when partons scattered at tt/2 were considered 
as moving in the medium with its rest system coinciding with the nuclei 
center of mass system. 

Two complications are to be considered. First, there are selected direc- 
tions of the current (z-axis) and its radiation. This is cured by introducing 
the permittivity tensor as 

D i {u,k) = e ij {u,k)E j {u,k). (87) 

Second, the impinging bunch content is similar to that of the target (plasma- 
plasma collisions!) and its permittivity must also be accounted for. In the 
projectile system it is the same as the permittivity of the target in its rest 



system, i.e. given by (66). Then one takes into account that the total induced 
current is the sum of currents induced in the target and in the projectile 
and Lorentz-transforms the projectile internal fields, polarization vector and 
currents to the target rest system where we consider the whole process. In 
this way the Lorentz transformation of the conductivity tensor <Jij(u>, k) (and, 
correspondingly, of e^u, k)) is found as it is done in [27)81 12091 I2TU1 12TT1 1212] . 



Using the current conservation and classical field equations (58), (59) with 
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g = one gets for non-zero components of the chromopermittivity tensor 



-yy 



-e. 



1 + + 

ur 7 



<Jih 



u 2 (u - k z )^' 
{u - k z ) 2 ^ 



(88) 



Here k^ and k z are the transverse and longitudinal components of k. We 
use the approximation of high energies (large 7 ^> 1 factor, i.e. v ~ 1). 
The terms depending on 7 are due to the impinging partons. They can be 
omitted everywhere except the terms which determine the Cherenkov gluon 
radiation at u — k z f» 0. 



The classical equations derived from (60), (61) and written in the mo- 



mentum space have solution if the following dispersion equation is valid 



det(u;, k) = \k 2 5ij — kikj 



0. 



(89) 



It is of the sixth order in momenta dimension. However, the sixth order 



terms cancel and (89) leads to two equations (of the second order): 



7 2 2 2 
k —00 — oo 



0. 



k- 



0. 



6/y 



(90) 
(91) 



They determine the internal modes of the medium and the bunch propagation 
through the medium, correspondingly. 



The equation (90 ) shows that the quark-gluon medium is unstable because 
there exists the branch with Im u > for modes k 2 < ojq. Thus the energy 
increase of the total cross sections responsible for positiveness of Re Fq(oj) 
at high energies is related to the instability of the quark-gluon medium. 



The equation (91) has solutions corresponding to Cherenkov gluons emit- 



ted by the impinging bunch and determined by the last term in (91). They 
can be found at u = k z + 5 (5 u>). For fc T = u one gets the solutions with 



Im 5\ 



2k z [2 1 (l + u 2 /k 2 )}^- 



(92) 
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For kf 7^ uq there is the solution with 



Im 5o 



Unk 



0*T 



kMkl-^ii + ujl/kDYn 



(93) 



It is well known (see |138j ) that the solutions of the disperion equation (89) 
determine the Green function of the system equations 



G(t,z) 



1 

2^ 



dk 



det(a;, k) 



exp(— iuit + ikz)du, 



(94) 



where the contour C(u) passes above all singularities in the integral. There- 



fore, the positive Im 5, in (92) and (93) correspond to the absolute instability 
of the system. Let us note that the instability at k T = u is stronger than 
at kx 7^ coo approximately by the factor 7 1 / 6 (this factor is about 4 times 



larger at LHC compared to RHIC). The instability exponent (92) decreases 
as 7 -1 / 3 and is about 16 times smaller at LHC compared to RHIC. It tends 
to zero asymptotically. 

Thus Cherenkov gluons are emitted with constant transverse momentum 
loq and their number, at high energies, is proportional to 



kn 



du 



where e(oo) is given by Eq. (66) with account of the threshold above which 
this equation is applicable. It differs from the traditional folklore of constant 
emission angle of Cherenkov radiation and the number of gluons oc du (or 
the total energy loss proportional to udu). This difference is easily explained 

const and the energy loss uxLj for e=const 

1 + LOI/UJ 2 . 



by Eqs (|70|), (|71|) which give cos 9 
and kx 



coo and dw/u energy loss for e 
If the data about Re Fq(uj) for hadronic processes are used in (|64|) and 



fitted by (66) then the value of uq depends on the density of scatterers N s 



and with above estimates can be of GeV order. 



4.5 Nonlinear effects and the color rainbow 



In the equations of in-medium QCD (|60j),(|61| considered in the subsection 



4.1 one assumed the chromopermittivity to be abelian, e° 



8' 



ab. 



The 



simplest generalization of the equations (60), (61) taking into account the 
nonabelian properties of the chromopermittivity e ab take into account the 
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tensor structure of e ab . Let us write the corresponding generalization of 
(60), (61) in the leading order in the coupling constant: 

.d 2 A b 



A A a -e 



ah 



dt 2 



e ab (A<5> b 



d 2 



e 6c I_<T) CN ) 
dt 2 ' 



-J 



(95) 



where (p a , j a ) are the components of the external color current j a defined in 
©■ 

From the physical meaning of the chromopermittivity e ab there follows 
that the most general expression for the chromopermittivity is a symmetric 
tensor e ab = e ba with the equal diagonal e' d ' and off-diagonal e < -°- ) components. 
To elucidate the physical content of the equations (95) let us consider the 
similarity transformation U diagonalizing the matrix e ab : 

= UeU' 1 . (96) 



In ( 96 ) the matrix e is the diagonal matrix with the eigenvalues of the chro- 
mopermittivity tensor {e'"'}, where a — 1, • • • , N 2 — 1, on the diagonal. It is 
easy to show that within the above-described assumptions on the structure 
of e ab one has: 



e( d ) + 2e(°), 
= - e< 



(97) 



The rotation Z7 allows to diagonalize (95) in terms of the modified fields 
= U^A^ and the modified current = U ab j b ^. 



AA"- f (fl) -A° 
ot 

9t 2 



9 2 r 



(98) 



From the above-described structure of the eigenvalues (97) there follows an 
existence of the two different solutions of the equations (98) corresponding 



7 In the present subsection we consider the simplest case of the real chromopermittivity. 
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to the eigenvalues (e*,e**): 



2g Q(vt - z - r ± y/e*v 2 - 1) 
e*" y/{vt - z) 2 - r 2 ± {e*v 2 -T) 
2g 0(vt — z — r±\/e**i 



$(DV-^-i( Tft ) = fj/ ^-'-'-Lv u tL , (99) 

e** ^/(!;t-z) 2 -ri(e**i; 2 -l) 

and the analogous solutions for A a . Returning to the original fields A a = 



~\ a 



U 1 A) one evidently gets the solution corresponding to the presence of the 



color Cherenkov rainbow formed by the Cherenkov radiation at the angles 
cos0* = ; 1 cos9** = . 1 (100) 

v ^ e (d) + 2e (o) vv / e (d) _ e (o) V ^ 



In the limit <^ the solution (99) turns into the standard solution 

To account for the nonlinear contributions to the in-medium Yang-Mills 
equations it is necessary to use the corresponding generalizations of (95). 
To demonstrate the possible practical consequences of such an approach in 
|207] one described the procedure of the model calculations of the nonlinear 
contributions to the gluon Cherenkov radiation, albeit with some simplified 
(and nonrealistic) assumptions concerning the color structure of the chro- 
mopermittivity. It was shown that in the case under consideration in the 
solutions there appear contributions of the form 

oc 



[(x-x ) 2 + (r,) 2 /4] 5/4 

and 



oc 



v^[(x-xo) 2 + (r t ) 2 /4] 3/4 ' 
which correspond to the color rainbow and differ (despite some similarity) 
from the lowest order solution ( 78 ) . 



4.6 Hydrodynamics (thermo dynamical and mechani- 
cal properties of QGP) 

The chromodynamical properties of QGP described above are its average 
characteristics measured by a test parton. Therefore they are more robust 



63 



than the thermodynamical and mechanical characteristics which change dur- 
ing the evolution of the quark-gluon medium. Spatial and temporal informa- 
tion is necessary to get the latter. The lattice calculations and ideas about 
CGC, Glasma, QGP provide some hints to transition from confined state to 
deconfined quarks and gluons. The inelastic collisions may lead to thermal- 
ization of the medium. Subsequently, it expands and hydrodynamics may be 
applied to treat this stage of evolution. It is described in many reviews (e.g., 
see |109[ IllOj ) and we give just a brief survey of it here. 

The medium is characterized by 6 independent variables. Those are the 
energy density e, pressure p, baryon number n B and 3 components of the 
velocity vector u^. The energy-momentum tensor and baryon number current 
are 

T^ix) = (e(x) + p(x))u»(x)u v (x) - p{x)g^, f B {x) = n B (x)u^(x). (101) 

The evolution of these variables is described by six equations of hydrody- 
namics. Those are five non-linear partial differential equations obtained from 
local conservation laws for energy, momentum and baryon number 

d lt T^(x) = (i/ = 0,..,3); dMx) = (102) 

and an equation of state relating p, e and tib- The latter is usually chosen in 
the closest possible form to the results of lattice QCD by normalization on 
the states below and above the critical temperature (i.e., hadrons to a quark- 
gluon medium). This is a rather arbitrary element of the whole approach. 
Besides, the solutions of the non-linear equations in (3+l)-dimensions ask for 
initial conditions to be defined and may be only obtained numerically with 
several external parameters. That is the origin of several conflicting results as 
was discussed above with references to [1251 EH EEl CEED HS21 US] . The 
CGC and Glasma approaches provide some guesses to the initial conditions 
for the evolution of a thermalized QGP. However, the evolution of the Glasma 
into a thermalized QGP is not yet understood. 

At least four parameters are necessary to fix the initial conditions and 
freeze-out algorithm which determines the transition from the hydrodynami- 
cal characteristics to hadronic stage. Those are the initial time r eq , energy (or 
entropy) density s eq , baryon number ns m and the freeze-out temperature (or 
decoupling energy density) e&, c . They are fixed by comparison of the theoret- 
ical results with the experimental data. Schematically, the correspondence 
of these parameters and experimental characteristics can be represented as 
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dN/dp T -T-e dec ; dN/dy- (rs) eq ; ■p/ir-n B ,eq/s eq \ (dN /dp T ) p / '(dN '/ 'dp T )^- 
r eq . Moreover, from the overlap geometry of the two colliding nuclei one can 
define the shape of the initial density profiles as well as make assumptions 
concerning the profiles of the initial longitudinal and transverse flows and 
the prescription for final hadronization. 

From solutions of the hydrodynamical equations one gets the transverse 
momentum spectra for various species of particles, the radial and elliptic 
flows, the shape of the interaction region (as revealed by Bose-Einstein corre- 
lations and Hanburry Brown-Twiss interferometry) . Their comparison with 
experimental data allows to determine the main thermodynamical, statis- 
tical and mechanical properties of the quark-gluon medium. Symbolically, 
their values found from 200 GeV data can be grouped as T eq « 360 MeV; 
T cr « 170 MeV; T dec « 120 MeV; r therm = r eq w 0.6 < 1 fm; r dec « 7 fm; 
e th « 25 GeV/fm 3 ; e cr w 1 GeV/fm 3 ; e dec « 0.075 GeV/fm 3 ; s eq w 110 fm" 3 ; 
rj/s ~ 0.1 (1/47T in AdS/CFT); n B < 0.5 fm -3 . These values show rapid 
thermalization, high average initial energy density and quite long "lifetime" 
of the quark-gluon plasma before hadronization at rather low energy density 
and temperature predicted by lattice QCD. Complete thermalization in the 
time less than 1 fm is required to obtain the measured value of elliptic flow 
and its centrality dependence which is very sensitive to any deviation from 
local thermal equilibrium of low px particles. Collective excitations, reso- 
nances and inelastic collisions keep the system in thermal equilibrium. The 
good agreement of the data with ideal fluid dynamics points to a very small 
viscosity of QGP. Other transport coefficients (shear, diffusion, heat conduc- 
tion) are not important if the microscale defined by rescattering is much less 
than the macroscale related to medium expansion. Strong non-perturbative 
interaction should be responsible for its behavior as an ideal liquid. Here- 
from comes the name of strong quark-gluon plasma (sQGP). In sQGP, there 
may exist, e.g., clusters |213j and colored bound states of massive quasi- 
particles [214] with heavier flavors. Resonance scattering on constituents of 
the quark-gluon medium can become important |215j . All that would give 
rise to its collective response with new scales discussed above, high pressure, 
large chromopermittivity and long-range correlations necessary to explain 
the enhancement of strange partons production by non-local processes. The 
long "lifetime" supports approximate description of energy losses within an 
"infinite" medium. 

Hydrodynamics is an actively developing field nowadays. Main charac- 
teristics of low pt particles have been described in this approach. At the 
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same time many factors have to be taken into account such as thermody- 
namics, models of collective flows, hadronization process, resonance decays, 
chemical composition, geometry of collisions etc. However a fully consis- 
tent hydrodynamical description has not yet been found. Some controversial 
conclusions appear from time to time about, e.g., the energy (from SPS 
to RHIC), rapidity and centrality dependence of elliptic flow, its absolute 
value, the transverse momentum dependence of various radii derived from 
HBT analysis, the chemical composition of some species (e.g., p/ir ratio) etc. 
One may hope that they will be resolved in the near future within the same 
set of adjusted parameters. That would provide deeper understanding of col- 
lective thermodynamical and mechanical properties of the bulk matter and 
its evolution. 



5 Some new possibilities at the LHC 

First LHC experiments on pp interactions at 2.36 and 7 TeV already showed 
that even the most general characteristics like inclusive pseudorapidity distri- 
butions in the central region \Arj\ < 2.4 for CMS (Compact Muon Solenoid) 
and | At/ | < 1 for ALICE (A Large Ion Collider Experiment) are not well 
enough described by theoretical predictions. Even though the qualitative 
features are reproduced, no Monte Carlo scheme (among the four tried) fits 
the results quantitatively. Surely, more surprises are waiting for us in corre- 



lation studies^ especially, in heavy-ion experiments. The energies of LHC 
open new possibilities for studies of the reviewed topics. The gluon-gluon 
processes will dominate. The formation length for radiation of gluons (and 
photons) becomes larger at higher energies. Thus high energy partons will 
feel the finite size of the interaction region. The spectra of radiated glu- 
ons and photons are predicted to become harder due to oscillations of the 
light-cone wave function [86| I216j . The emission from heavy quarks should 
be enhanced by the finite-size effects relative to light quarks. In general, 
the relative roles of surface effects, initial high virtualities of created high-p-^ 
partons, dead-cone effect in radiation by heavy quarks will be elucidated. 
Coherent final state interactions must become stronger. 

The smaller values of chromopermittivity for very high energy partons 



18 The CMS Collaboration reported |218| an observation (!) of a ridge with the narrow 
distribution in relative azimuthal angle and broad distribution in pseudorapidity difference 
in pp collisions at the energy of 7 TeV. 
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discussed above indicates on weaker coupling and lower parton density of 
the medium at LHC compared with RHIC. Therefore, the role of fluctuations 
induced by instabilities will be enhanced at LHC |217j . 

Both trigger and non-trigger experiments are important in AA-collisions. 
Following RHIC in trigger experiments, the partons (jets) at large angles with 
much higher energies will be available. Correspondingly, this allows to study 
jet quenching, ridge and double-humped events in wider energy intervals. In 
particular, the energy dependence of the chromopermittivity may be found 
by studying humps behavior in wider energy intervals. The humps shifts due 
to the wake effect in semi-central collisions and due to Lorentz transformation 
of the chromopermittivity tensor in a wide energy interval are of a special 
interest. 

In non-trigger experiments studies of the response of the quark-gluon 
medium to the extremely fast moving partons will be of primary importance. 
It would be quite instructive if the events similar to the cosmic ray one are 
observed. The finite-size effects are especially important for forward-moving 
partons with large x. 

The search for forward subjets with a fixed mass (due to the transition 



radiation - Eq. (86)) or fixed transverse momentum (due to high energy 



Cherenkov gluons - Eq. (91)) might be fruitful. 

Effects due to Cherenkov gluons can be used to study the wide region of 
(x, Q 2 )-plane available at the LHC. The different values of parton densities 
would correspond to different subregions on this plane. 

The nonlinear effects described in section 4.5 may become more pro- 
nounced albeit still difficult to separate from main effect due to the possible 
dispersive behavior of the chromopermittivity. 

One can imagine more exotic experiments. For example, let triggers reg- 
ister a single electron (positron) which is produced in a process of the convo- 
lution of the quark- ant iquark pair inverse to those studied at e + e~-colliders. 



The dielectric permittivity can be measured and the formulae (63), (66) com- 
pared with the data. Even though the cross section is very small, it might 
be measurable in the high luminosity regime of LHC. 

6 Conclusions 

The experimental studies at SPS and, especially, the numerous and detailed 
data about heavy ion collisions at RHIC provided us with invaluable infor- 
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mation about the completely new field - the properties of the quark-gluon 
medium. They showed that AA events can not be described as an addi- 
tive collection of pp collisions. The heavy ion dynamics must take into ac- 
count the collective behavior of the medium. That is clearly demonstrated 
by anisotropic flows, jet quenching, such peculiar correlations as ridge and 
double-humped events etc. 

Theoretical understanding of medium evolution asked for application 
of QCD at full strength to speculate about such stages as CGC, Glasma, 
thermalization, QGP, hadronization. The specific configuration of fields in 
Glasma and corollaries of deconfinement at QGP stage need further studies. 
In its turn, experimental results required to revive the theoretical methods 
widely applied in the condensed matter physics. The modification of en- 
ergy losses of partons in the medium because of the formation length effects 
and its collective response (the chromopermittivity) due to polarization are 
described in this manner. Hydrodynamics was extensively exploited for de- 
scription of collective behavior of this medium. 

Both experimental and theoretical progress in this field are extremely 
fruitful and promising. Surely, even more interesting results will be obtained 
in further searches at RHIC. Exciting perspectives are opened with the advent 
of LHC in the new energy region. 
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Figure captions. 

Fig. 1. QCD cascade in the medium. 

Fig. 2. Distribution in rapidity P(y) of final prehadrons: 1) L — Ofm, 
full angular ordering, red, full line; 2) L = 0.5 fm, partial angular ordering, 
green, dashed line; 3) L — 5fm, partial angular ordering, blue, dotted line 
[97j. 

Fig. 3. Associated azimuthal correlations at STAR: experiment- circles, 
theory- triangles. 
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Fig. 4. Modified spectrum of dileptons in semi-central collisions In(158 
A GeV)-In for NA60 (points) compared to the p-meson peak in the medium 
with additional Cherenkov contribution (dashed line). 

. 5. Distribution in the number of produced particles at different dis- 
tances from the collision axis r in the stratospheric event at 10 16 eV [T37] 
has two pronounced peaks. Correspondingly, the distribution in pseudora- 
pidity also has two such peaks. 
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